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Abstract 

An  understanding  of  spontaneous  emission  processes  within  microcavities  is 
crucial  in  addressing  the  need  to  make  tomorrow’s  microlasers  more  efficient.  One 
approach  to  improving  the  device  efficiency  is  to  reduce  the  threshold  input  current 
at  which  lasing  begins  to  occur.  It  has  been  suggested  that  the  threshold  in  a  mi¬ 
crocavity  laser  can  be  decreased  by  increasing  the  fraction  of  spontaneous  emissions 
into  the  lasing  mode.  This  can  be  accomplished  by  preferentially  coupling  the  gain 
medium  of  the  laser  to  the  electromagnetic  cavity  mode  of  interest.  It  therefore 
becomes  necessary  to  understand  the  mechanism  by  which  this  coupling  takes  place. 
This  research  develops  a  fully  quantum  mechanical  description  of  the  interaction  be¬ 
tween  a  gain  medium  modeled  as  a  two  level  atom  and  a  multimode  electromagnetic 
field  in  a  microcavity.  Atomic  transition  probabilities  are  computed  for  systems  in 
which  the  atom  couples  through  a  single  photon  process  to  electromagnetic  cavity 
modes  which  range  in  number  from  two  to  2000.  Calculations  performed  for  cavities 
with  widely  spaced  modes  demonstrate  that  atoms  exhibit  Jaynes- Cummings  behav¬ 
ior  when  closely  tuned  to  one  mode.  Detuning  of  the  atom  from  the  mode  inhibits 
the  exchange  of  energy,  while  increasing  the  strength  of  the  coupling  to  the  mode 
amplifies  this  exchange.  Two  level  systems  strongly  coupled  to  many  closely  spaced 
modes  exhibit  spontaneous  emission  rates  characteristic  of  an  atom  in  free  space. 


XIV 


SPONTANEOUS  EMISSION  IN  MICROCAVITY  LASERS 


I.  Introduction 

1.1  Motivation 

As  stated  in  the  New  World  Vistas,  one  of  the  goals  of  the  Air  Force  laser 
programs  heading  into  the  twenty-first  century  is  to  develop  very  compact,  highly 
efficient  semiconductor  laser  devices  [24].  The  breadth  of  military  applications  for 
miniaturized  laser  devices  continues  to  grow  as  portability  and  cost-efficiency  of 
mass-manufcicturing  become  increasingly  important  features  across  a  variety  of  disci¬ 
plines:  optical  control  circuitry,  information  storage  and  retrieval,  optical  encryption 
and  decryption  technologies — just  to  name  a  few.  In  the  past,  the  applications  for 
semiconductor  lasers  were  satisfied  almost  exclusively  by  traditional  edge-emitting 
devices.  These  conventional  “horizontal”  lasers  range  in  length  from  roughly  200  to 
450  /am,  though  shorter  devices  have  been  fabricated  [23].  The  design  is  straight¬ 
forward:  an  active  region  is  formed  by  layering  appropriate  junction  materials  that 
also  act  as  the  optical  waveguide  in  operation.  The  ends  of  the  device  are  cleaved 
or  etched  to  form  mirror  surfaces  which  are  orthogonal  to  the  junction  (see  Figure 
1.1a).  In  this  way,  the  cavity  of  the  resonator  is  along  the  interface  of  the  active 
region  (horizontal).  In  these  Fabry  Perot  resonators,  relatively  long  cavity  lengths 
are  required  to  reduce  the  effects  of  loss  at  the  mirrors.  However,  the  resonant 
mode  spacing,  or  Free  Spectral  Range,  is  inversely  proportional  to  the  cavity  length, 
which  in  turn  means  long  cavities  support  nearly  a  continuum  of  resonant  modes 
which  compete  for  the  energy  released  by  the  active  region  at  the  junction.  What  is 
needed  is  a  means  of  taking  the  non-lasing  modes  out  of  the  competition. 
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Active  Region 


Mirror  Surface 


b) 

Figure  1.1  Diagrams  for  (a)  conventional  edge-emitting  diode  laser  and  (b)  VCSEL 
diode  laser  (After  Ref.  [23]). 


1.2  Emergence  of  Vertical  Cavity  Surface  Emitting  Lasers 

Within  the  past  decade,  research  in  the  area  of  alternate  semiconductor  laser 
designs  has  produced  the  Vertical  Cavity  Surface  Emitting  Laser  (VCSEL).  In  con¬ 
trast  to  the  edge-emitting  approach,  the  VCSEL  design  utilizes  a  layering  technique 
to  form  Distributed  Bragg  Reflectors  (DBR)  at  either  end  of  the  cavity  such  that 
the  mirrors  of  this  geometry  are  parallel  to  the  active  region  in  the  middle  (Figure 
1.1  (b)).  These  mirrors  axe  formed  by  layering  pairs  of  alternately  low  and  high 
refractive  index  materials  as  shown  in  Figure  1.2. 

VCSELs  are  fabricated  using  a  bottom-up  layering  process  using  techniques 
such  as  Molecular  Beam  Epitaxy  or  Metal- Organic  Chemical  Vapor  Deposition  to 
grow  the  entire  device  on  the  substrate  from  the  first  DBR  to  the  active  region  to  the 
second  DBR.  Although  tolerances  on  the  layering  process  are  extremely  tight,  this 
design  eliminates  the  often  complex  cleaving  or  etching  steps  necessary  to  produce 
the  mirrors  in  the  conventional  laser  diodes.  In  the  final  product,  VCSEL  cavities 
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measure  in  the  micron  to  submicron  range,  roughly  two  orders  of  magnitude  shorter 
than  edge-emitting  diodes.  The  advantages  of  this  technology  over  the  edge-emitting 
device  are  well  documented  and  include  single  longitudinal  mode  operation  (leading 
to  higher  output  in  the  desired  mode),  smaller  threshold  currents,  and  the  potential 
for  application  in  large  two-dimensional  laser  arrays  [3,23].  In  what  way  can  these 
devices  be  made  more  efficient?  It  has  been  suggested  by  Vurgaftman  and  Singh 
that  the  threshold  current  for  microcavity  lasers  can  be  reduced  even  further  through 
reductions  in  cavity  dimensions  [25].  It  is  instructive  at  this  point  to  momentarily 
depart  from  the  present  discussion  to  address  the  basic  principles  of  operation  for 
laser  diodes  then  return  to  the  implications  of  cavity  length  reductions. 

1.2.1  Overview  of  the  VCSEL.  At  the  risk  of  oversimplifying,  the  VCSEL 
can  be  described  in  terms  of  three  basic  components: 

•  an  active  gain  medium; 

•  the  pump  energy,  typically  in  the  form  of  an  applied  electrical  current  or  optical 
pumping  that  excites  the  active  medium; 

•  an  optical  cavity  formed  by  the  DBR  mirrors  that  encloses  the  active  medium 
and  allows  for  oscillation  and  subsequent  amplification  of  emitted  radiation; 

It  is  the  interactions  between  the  active  medium  and  the  electromagnetic  field  in  the 
cavity  with  which  this  thesis  is  chiefly  concerned. 

1.2.2  Emission  Processes  Within  The  Active  Medium.  In  very  elementary 
terms,  the  active  medium  consists  of  the  region  surrounding  an  interface  between 
semiconductor  materials  and  is  chajacterized  by  an  energy  bandgap  structure.  A 
simple  structure  is  shown  in  Figure  1.3.  In  this  bandgap  structure,  charge  carriers  are 
excited  from  a  lower  energy  bound-state  valence  band  into  a  higher  energy  conduction 
band  by  the  pump  energy.  Once  excited,  the  carriers  may  transition  back  to  the 
valence  band,  either  directly  or  by  way  of  intermediate  transitions.  The  desirable 
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Transmitted 
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Incident 

Figure  1.2  Multilayer  reflections  of  the  electric  field  within  the  DBR. 

transitions  are  those  for  which  the  de-excitation  is  accompanied  by  the  emission 
of  photons  of  the  intended  lasing  wavelength;  other  means  of  decay  detract  from 
the  overall  efliciency  of  the  device  and  require  higher  pumping  levels  to  overcome 
this  effect.  When  a  radiative  decay  is  ‘uninduced’,  or  independent  of  the  photon 
population  in  the  cavity,  the  process  is  referred  to  as  spontaneous  emission.  When 
a  photon  is  emitted  into  the  cavity,  there  is  also  a  probability  for  it  to  interact 
with  other  carriers  in  the  region  in  one  of  two  ways:  through  absorption.,  where  the 
photon  is  absorbed  to  excite  the  carrier,  and  stimulated  emission,  where  the  passage 
of  a  photon  induces  the  excited  carrier  to  de-excite  and  emit  a  photon  of  matching 
wavelength,  polarization,  phase,  and  direction.  If  a  population  inversion  is  achieved 
and  maintained  such  that  the  number  of  conduction  band  carriers  exceeds  that  in 
the  valence  band,  then  the  subsequent  avalanche  of  spontaneous  emissions  create  a 
substantial  flux  of  photons  in  the  cavity.  A  fraction  of  these  photons  are  emitted  into 
a  resonant  mode  of  the  cavity  and  remain  in  the  cavity  temporarily  until  absorbed 
or  lost  as  output.  While  in  the  cavity  the  photons  can  initiate  stimulated  emission. 
When  the  gain  in  stimulated  emission  exceeds  loss  to  output  or  absorption,  then 
lasing  threshold  is  achieved. 


E 


Figure  1.3  Simple  energy  bandgap  structure. 

1.3  Modified  Spontaneous  Emission 

We  now  briefly  return  to  the  idea  of  modifying  the  probability  of  spontaneous 
emission  into  the  lasing  mode.  In  order  to  improve  the  efficiency  of  the  laser,  it  is 
useful  to  reduce  the  power  required  to  bring  the  laser  to  threshold.  One  method  of 
accomplishing  this  is  to  preferentially  couple  the  gain  medium  to  the  lasing  mode  and 
decrease  the  coupling  to  other  modes.  It  has  been  suggested  that  the  degree  to  which 
the  medium  couples  to  a  mode  can  be  enhanced  by  tuning  the  transition  frequency 
in  the  medium  to  the  mode  frequency  [28].  The  Jaynes-Cummings  Model,  proposed 
in  1963,  demonstrates  that  a  two  level  system  interacting  with  a  single  mode  with 
which  it  is  in  resonance  will  exchange  energy  completely  with  that  mode  [13]. 

It  is  in  this  context  that  the  research  described  in  this  thesis  is  presented.  A 
theoretical  model  is  presented  which  treats  the  electromagnetic  field,  gain  medium, 
and  their  coupling  within  a  microcavity  laser  quantum  mechanically.  This  model  is 
an  extension  of  the  Jaynes-Cummings  theory  of  the  interaction  between  an  atomic 
spontaneous  emission  source  and  the  vacuum  field  mentioned  above. 
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Figure  1.4  (a)  Spontaneous  Emission,  (b)  Absorption,  and  (c)  Stimulated  Emis¬ 

sion. 

1.4  Jaynes-Cummings  Model 

The  Jaynes-Cummings  Model  (JCM)  as  presented  by  Meystre  and  Sargent  [17] 
models  the  interaction  between  an  atom  represented  by  a  two  level  system  and  a 
single  electromagnetic  mode  in  an  enclosing  optical  cavity.  The  model  assumes  a 
single  photon  mechanism  in  which  a  single  quantum  of  energy  is  either  absorbed 
from  the  field  by  an  atom  in  the  ground  state  or  emitted  into  the  field  mode  by 
an  excited  atom  The  time  evolution  of  the  system  is  then  given  in  terms  of  the 
evolution  of  the  probability  for  the  atom  to  spontaneously  emit  into  the  field  and 
decay  to  the  ground  state.  The  JCM  predicts  a  sinusoidal  interchange  of  probability 
between  the  excited  atom  {P  [i))  and  excited  field  (P  {2))  states  [17]: 

-  S^  +  4g^  +  +  (1-1) 

where  p  is  a  constant  that  describes  how  strongly  the  atom  is  coupled  to  the  field 
mode  and  6  is  the  detuning  of  the  atomic  transition  frequency  from  the  field  mode. 

^The  definition  of  the  term  ‘photon’  varies  widely  according  to  the  context  in  which  it  is  [)re- 
sented.  In  this  thesis,  ‘photon’  is  used  to  indicate  a  single  quantum  of  electromagnetic  energy  in  a 
particular  mode  of  the  cavity.  An  interesting  discussion  of  the  definitions  of  the  term  is  found  in 
Ref.  [17]  in  Section  12-1. 
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The  resulting  behavior  of  these  probabilities  is  known  as  Rabi  flopping  after  LI.  Rabi 
(1936),  who  explored  the  behavior  of  spin-|  magnetic  dipoles  in  nuclear  magnetic 
resonance  [17].  The  frequency  Tl  =  +  4^2  jg  the  Rabi  frequency  for  systems 

with  no  photons  in  the  cavity  when  the  atom  is  excited.^  The  probability  exchange 
described  by  Equations  1.2  is  illustrated  in  Figures  1.5a-d  and  1.6a-d. 

It  is  clear  that  as  the  coupling  grows  stronger,  the  frequency  of  the  flopping 
increases,  indicating  a  more  rapid  transition  back  and  forth  between  the  upper  and 
lower  atomic  states.  When  the  atom  and  field  are  perfectly  tuned,  complete  exchange 
occurs.  An  increase  in  8  increases  the  frequency,  but  also  reduces  the  amplitude 
of  the  exchange  so  that  the  energy  in  the  system  increasingly  remains  in  the  atom. 
However,  for  g  ^  5,  the  exchange  amplitude  for  a  detuned  atom  approaches  1.0.  The 
JCM  demonstrates  that  the  degree  of  exchange  between  the  atom  and  field  mode  is 
suppressed  through  detuning  and  enhanced  for  increasing  coupling  strength. 

As  mentioned  previously,  the  mode  spacing,  or  Free  Spectral  Range  (FSR),  in 
the  cavity  is  inversely  proportional  to  the  cavity  length.  For  macroscopic  cavities 
such  as  those  used  in  typical  gas  discharge  lasers,  the  mode  spacing  is  very  small. 
An  atom  near  resonance  with  one  mode  is  likely  as  well  to  be  near  resonance  with 
adjacent  modes.  Conversely,  as  the  cavity  length  is  reduced,  the  mode  spacing 
increases.  In  the  limit  of  microcavity  structures  (roughly  on  the  order  of  tens  of 
microns  or  smaller),  the  mode  spacing  becomes  very  large  (consider  that  the  FSR 
for  a  10  fim  cavity  is  on  the  order  of  10^^  Hz  compared  to  10®  Hz  for  a  1  m  cavity). 
An  atom  neax  resonance  with  a  given  mode  in  a  microcavity  is  much  more  detuned 
from  the  adjacent  modes.  Therefore,  in  the  context  of  the  microcavity,  the  single 
mode  description  contained  in  the  JCM  becomes  applicable. 

^The  more  general  form  of  the  Rabi  frequency  is  7?  =  +  4  (ra  +  1)  where  n  is  the  number 

of  photons  in  the  cavity  when  the  atom  is  excited. 
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Probability  Probability 


Figure  1.6  Increasing  the  coupling  between  the  atom  and  the  single  field  mode 
increases  the  fiopping  frequency  but  does  not  affect  the  amplitude  in 
the  JCM.  5  =  0  (a)  Uncoupled  (b)  go  (c)  4  go  (d)  8  go 
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1.5  Scope 

A  truly  quantum  mechanical  description  of  the  microlaser  system  as  a  whole 
requires  a  quantized  description  of  every  facet  of  the  problem.  The  individual  atoms 
in  the  gain  medium  are  in  and  of  themselves  complex  systems.  Because  they’re 
not  isolated  from  one  another,  interactions  occur  which  must  be  considered  as  well. 
Similarly,  the  atoms  in  the  cavity  walls  must  interact  among  themselves  and  with 
the  gain  medium,  and  the  whole  collection  of  material  interacts  with  the  field  in  the 
cavity.  Broadening  the  field  of  view  brings  external  influences  into  the  picture.  Each 
component  of  this  system  can  taJce  on  a  quantum  mechanical  description,  and  the 
resulting  problem  grows  well  beyond  the  number-crunching  capabilities  of  even  the 
best  algorithms  and  high-performance  computers.  In  order  to  make  this  problem 
tractable,  a  number  of  simplifications  wUl  be  made: 

•  As  in  the  JCM,  the  gain  medium  will  be  modeled  by  a  single  two-level  atom. 
In  etfect,  this  corresponds  to  taking  a  vertical  slice  through  the  band  structure 
shown  in  Figure  1.3.  In  principle,  the  medium  is  best  modeled  by  a  continuum 
of  these  slices,  but  the  two  level  representation  will  have  to  suffice  in  this 
simplified  model. 

•  A  quantum  mechanical  description  of  the  cavity  wall  material  is  not  attempted 
and  the  walls  are  assumed  to  be  perfectly  electrically  conducting  (PEC).  As  a 
result,  the  only  interactions  axe  between  the  atom  and  the  field,  and  there  is 
no  loss  to  the  cavity  walls. 

•  Because  the  wavelengths  under  consideration  fall  in  the  optical  regime  and  are 
therefore  roughly  a  thousand  times  larger  than  the  dimensions  of  the  atom  in 
the  cavity,  the  atom/field  interaction  will  be  approximated  as  the  interaction 
of  an  electric  point  dipole  in  an  electromagnetic  field.  Magnetic  dipole,  electric 
quadrupole,  and  higher  multipole  terms  will  be  neglected. 
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•  The  two-level  system  is  assumed  to  be  stationary  with  regard  to  position  or 
moving  slowly  enough  within  the  cavity  that  variations  in  positions  are  small 
enough  to  be  neglected  on  the  timescale  of  the  interactions.  This  carries  the 
added  benefit  of  allowing  one  to  eliminate  velocity  dependent  Doppler  shifts 
from  the  picture. 

•  A  one  dimensional  treatment  of  the  electromagnetic  field  was  adopted  to  sim¬ 
plify  calculations.  The  results  may  be  extrapolated  to  the  three  dimensional 
cubic  cavity  without  significant  extra  work,  where  the  spatial  dependence  of 
the  electric  and  magnetic  fields  is  slightly  a  slightly  more  complex  product  of 
sines  and  cosines  in  x,  y,  and  z.  In  addition,  the  mode  frequencies  lose  their 
harmonic  spacing  in  the  translation  to  the  three  dimensional  cavity.  ^ 

•  Cavity  energies  are  low  enough  to  warrant  a  strictly  non-relativistic  treatment. 

These  assumptions  and  constraints  form  the  framework  within  which  the  following 
approach  is  pursued. 

1.6  Methodology 

As  mentioned  earher,  the  model  to  be  developed  is  an  extension  to  the  JCM 
outlined  in  Section  1.4.  Where  the  JCM  coupled  the  atom  to  a  single  mode  of  the 
electromagnetic  field  within  the  cavity,  here  the  freedom  to  select  multiple  modes  is 
added.  In  the  model  developed  in  this  thesis,  the  atomic  system  is  quantized  in  a 
fashion  completely  analogous  to  the  spin-|  particle  in  a  magnetic  field.  Transitions 
between  the  two  atomic  levels  are  accomplished  by  means  of  operation  on  the  state 
vector  with  the  atomic  raising  and  lowering  operators  (see  Chapter  2).  Quantization 
of  the  electromagnetic  field  is  accomplished  by  first  defining  the  mode  structure  by 
solving  the  classical  cavity  problem  for  a  PEC  cavity,  then  translating  to  a  quantum 

^The  1-D  frequencies  vary  linearly  with  mode  number  m  whereas  in  a  3-D  cubic  cavity  the 
dependence  is  (n^  -f  -1- 

^For  an  excellent  discussion  of  spin-j  systems,  see  Cohen-Tannoudji  et  al,  Ch.  IV  (Ref.  [8]) 
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mechanical  picture  by  associating  with  each  mode  the  characteristics  of  a  quantized 
harmonic  oscillator.  A  dipole  interaction  is  assumed,  and  coupling  coefficients  are 
defined  which  determine  the  strength  of  interaction  between  the  atom  and  each  of  the 
cavity  modes.  Together,  the  atomic,  field,  and  coupling  terms  yield  the  full  Hamil¬ 
tonian  of  the  system.  The  full  Hamiltonian  is  diagonalized  and  the  time  evolution 
of  the  eigenstates  is  oscillatory  with  characteristic  frequencies  given  by  the  energy 
eigenvalues.  Finally,  the  transformation  matrix  that  diagonalizes  the  Hamiltonian  is 
used  to  transform  back  into  the  uncoupled  basis  to  determine  as  a  function  of  time 
the  probability  that  the  atom  has  undergone  spontaneous  emission. 

The  results  obtained  in  Chapter  6  were  obtained  by  computing  the  matrix 
elements  of  the  full  Hamiltonian  for  the  system  and  diagonalizing  the  resulting  ma¬ 
trix  using  FORTRAN  90  and  a  public  domain  Linear  Algebra  Package  (LAPACK) 
diagonalizing  routine.  The  resulting  eigenvectors  and  eigenvalues  were  then  used  to 
propagate  the  system  in  time  giving  the  evolution  of  the  probabilities  to  be  in  each 
of  the  uncoupled  states. 

1.7  Notation 

Where  applicable  I’ve  adopted  the  Dirac  Bra-Ket  notation.  State  vectors  are 
represented  by  |  9?)  €  V,  where  V  is  a  vector  space.  Additionally,  \cip)  =  c  |  9?) 
with  c  in  general  a  complex  constant.  The  scalar  product  is  defined  as  {ipi\(p2)  =  c 
with  the  corresponding  complex  conjugate  (^2|v^i)  =  c*.  Abstract  operators  are 
designated  by  a  “hat”  above  the  operator:  A,  for  example.  The  hat  is  dropped  when 
a  representation  is  taken,  as  with  (a^  |  A  |  ccy),  where  aij  are  members  of  the  basis 
in  which  A  is  represented.  Note  however,  that  x,  y,  and  z  are  used  as  unit  vectors 
in  the  classical  field  discussions  to  denote  the  direction  of  a  vector  quantity  with 
respect  to  a  set  of  x,  y,  and  z  coordinate  axes,  such  as  the  classical  electric  field 
vector  IEj  =  Eqx  directed  along  the  x  direction.  Finally,  the  matrix  elements  of  the 
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operator  A  are  Aij  =  {ui\  A  \uj)  such  that 
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1.8  Overview 

The  remainder  of  this  thesis  is  devoted  to  presenting  the  theory,  computational 
results,  and  implications  of  this  quantum  mechanical  approach  to  the  microcavity 
laser.  To  this  end,  a  chapter  has  been  devoted  to  each  of  these  three  aspects  of  the 
problem  as  outlined  in  the  following  paragraphs. 

The  focus  of  Chapter  2  is  to  present  the  theory  underlying  the  quantum  me¬ 
chanical  treatment  of  the  atom.  The  analogy  is  made  between  the  two  level  system 
and  a  spin-|  particle  by  means  of  the  atomic  raising  and  lowering  (or  spin- flip)  op¬ 
erators.  The  resulting  atomic  Hamiltonian  is  then  represented  in  the  |  ±)  eigenbases 
in  matrix  format. 

In  Chapter  3,  expressions  for  the  electric  and  magnetic  fields  are  obtained 
for  the  classical  PEC  cavity.  The  amplitudes  of  these  quantities  are  normalized 
by  considering  the  total  energy  in  the  cavity,  and  the  translation  to  a  quantum 
mechanical  operator  form  follows.  As  in  Chapter  2,  a  matrix  form  is  obtained  by 
representing  the  Hamiltonian  in  the  eigenbasis,  in  this  case  where  the  {1,} 

notation  indicates  that  the  state  has  a  single  photon  in  the  mode  and  none 
in  any  of  the  others.  The  uncoupled  Hamiltonian  is  then  obtained  from  the  sum  of 
the  atomic  and  field  Hamiltonians. 

Chapter  4  develops  the  dipole  interaction  between  the  atom  and  field.  The 
interaction  Hamiltonian  matrix  is  obtained  by  representing  the  Hamiltonian  in  the 
uncoupled  atom/field  basis  {  |  -H  {0}),  |  —  {!,})  }. 
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The  atom,  field,  and  interaction  Hamiltonians  are  brought  together  in  Chapter 
5.  The  chapter  discusses  how  the  eigenvectors  of  the  diagonalized  Hamiltonian  are 
used  to  translate  between  the  eigenbasis  and  the  uncoupled  basis  mentioned  above. 
The  time  evolution  operator  is  introduced  here,  and  the  steps  required  to  compute 
the  time  evolution  of  the  bare  state  probabilities  are  outlined. 

Chapter  6  contains  the  results  of  the  calculations  based  on  the  theory  presented 
in  Chapters  2  through  5.  A  number  of  parameters,  including  cavity  size,  atomic 
transition  frequency,  and  the  coupling  factor,  are  varied  and  the  eflFects  of  these 
changes  on  the  evolution  of  the  system  are  extracted.  Two  regimes  are  of  particular 
interest: 

•  atoms  interacting  with  a  field  of  widely  spaced  modes  under  weak  coupling 

•  atoms  strongly  coupled  to  a  large  number  of  tightly  spaced  modes 

In  the  limit  of  weak  coupling  and  small  cavities,  the  results  obtained  reproduce 
those  of  the  JCM  oscillations  with  predicted  Rabi  flopping  frequencies.  Likewise,  for 
stronger  coupling  and  large  cavities  where  the  atom  easily  accesses  a  large  number  of 
modes,  the  atom  decays  exponentially  from  the  excited  state  with  a  rate  comparable 
to  an  atom  in  free  space.  Chapter  6  concludes  with  a  brief  discussion  of  several 
unanticipated  results. 

The  final  chapter  in  this  thesis.  Chapter  7,  presents  a  summary  of  the  results 
and  potential  implications  for  microlciser  performance.  Recommendations  on  future 
studies  in  this  area  are  provided  at  the  closing. 

Appendix  A  will  describe  the  steps  taken  to  verify  the  validity  of  the  FOR¬ 
TRAN  code  written  for  this  problem,  while  the  code  itself  is  presented  as  Appendix 
B.  Appendix  C  contains  several  plots  demonstrating  an  anomalous  suppression  of 
atomic  de-excitation  for  moderate  coupling  strengths.  Finally,  Appendix  D  lists  sev¬ 
eral  useful  conversion  factors  from  the  atomic  unit  system  to  the  MKS  and  Gaussian 
unit  systems. 


1-14 


11.  The  Atomic  Hamiltonian,  Ha 

This  chapter  develops  the  Hamiltonian  describing  the  atomic  system.  Recall  from 
Section  1.5  that  the  simplification  is  made  to  represent  the  atom  as  a  two- level 
system.  While  this  significantly  simplifies  this  aspect  of  the  model,  it  is  by  no  means 
intended  as  a  realistic  model  of  the  material  within  a  VCSEL  or  other  microcavity 
laser. 

2.1  Two  Level  System 

Consider  the  two  level  system  shown  in  Figure  2.1.  The  upper  and  lower  states 
of  the  system  are  written  |-|-)  and  |— )  and  are  eigenstates  of  the  atomic  Hamiltonian, 
Ha: 

Ha  |±)  =  |±) 

with  energy  eigenvalues  hu)±.  Note  that  the  operator  Ha  can  also  be  written 


Ha  =  1  Ha  1 


where  1  is  the  unit  operator.  The  closure  relation  of  quantum  mechanics  becomes 
useful  at  this  point: 

E  i*X'i  =  1 

i 

if  the  |i)  form  a  complete  orthonormal  basis.  The  previous  expression  for  Ha 
becomes 


Ha  =  E  l‘H‘l  I E  WO'I 

«  3 

=  'Z\i)'£{i\HM{j\ 

*  3 

• 

i  j 
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Figure  2.1  Renormalization  of  the  atomic  energy. 


Recall  that  (f|j)  =  6ij  for  orthonormal  bases.  Therefore, 

Ha  = 

a 

=  ^ha;i|i)(i|.  (2.1) 

i 

For  the  convenience  of  subsequent  discussion,  suppose  the  zero-energy  reference  is 
renormalized  to  a  value  precisely  halfway  between  the  energies  of  the  two  states  as 
in  Figure  2.1.  If  the  energy  of  transition  between  the  two  levels  is 


hu)j^.  —  hw-  =  —  U-)  =  hu  , 

then  renormalization  of  the  energies  presents  a  new  set  of  measured  values 

/  I  rr  I  \  hu;  ,  ,  .  hU> 

(+|ff^l+>  =  =  T 


The  operator  Ha  can  finally  be  re-expressed  (with  the  help  of  Equation  2.1)  as 


=  yd+X+l-l-X-l). 


(2.2) 


2.2  Matrix  Representation  of  Ha 

At  this  point  it  is  worthwhile  to  discuss  the  matrix  representation  of  the  atomic 
Hamiltonian.  Consider  the  expansion  of  in  the  |±)  basis: 


l+X+l 

(-1 


1+)  |-> 

(  (+l+>(+l+>  (+l+>{+|->  \ 
V  (-i+>(+i+>  (-I+X+I-) ) 


(2.3) 


Those  elements  containing  the  scalar  product  (±|t)  vanish  by  orthogonality,  while 
the  remaining  element  is  1,  leaving 


l+X+l « 


This  matrix  can  be  written  as  the  sum  of  two  others: 


l+)(+l  ^ 


(2.4) 


The  first  matrix  in  the  brax:kets  is  the  identity  matrix,  1,  while  the  second  is  the 
matrix  representation  of  the  Pauli  spin  operator.  Performing  similar  operations 
for  j— )(— I  leads  to  the  correspondences 


l+X+l 

l-x-l 


1  +  (Tz 


1-0-^ 


(2.5) 

(2.6) 
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When  the  results  of  Equation  2.6  are  substituted  into  Equation  2.2,  an  expression 
for  Ha  is  obtained  in  terms  of 


1  +  cr 2 
2 


1^) 


And  in  operator  form, 


Ha 


(2.7) 
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III.  The  Multimode  Electromagnetic  Field  Hamiltonian,  Hp 

Attention  turns  now  to  developing  the  Hamiltonian  for  the  electromagnetic  field  in 
the  cavity.  Following  the  method  of  Chapter  2,  the  ultimate  goal  here  is  an  expression 
for  Hp,  the  field  Hamiltonian,  in  the  bare  state  representation.  Initially,  the  classical 
modes  of  the  cavity  are  derived  in  Gaussian  units  with  a  transition  to  atomic  units  for 
the  quantum  mechanical  expressions.  Throughout  the  classical  cavity  development, 
the  quantities  cq  and  yUo,  which  are  both  equal  to  1  in  these  units,  are  retained  as 
place  holders  to  facilitate  future  studies  that  investigate  dispersive  cavities. 

3.1  Classical  Electromagnetic  Cavity  Modes. 

The  classical  cavity  discussion  begins  with  Maxwell’s  equations  for  free  space 
in  Gaussian  units, 


V,  T.  1 

V  X  E+ =  0 
c  at 


a) 


c  dt 


b) 


VE  =  0 


c) 


VB  =  0 


d) 


The  curl  of  the  first  equation  yields 

1  a 

VxVxE  +  -— (VxB)  =  0. 
c  at 


(3.1) 


(3.2) 
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Figure  3.1  One  dimensional  PEC  microcavity. 


Substituting  from  Equation  3.1b  and  expanding  the  repeated  curl  above  produces 
the  wave  equation  in  E: 


V^E- 


c2  df^ 


(3.3) 


This  equation  can  now  be  subjected  to  the  boundary  conditions  for  the  PEC  cavity. 
Choosing  the  z-axis  as  the  direction  along  the  cavity  length  and  an  E  polarization 
in  the  x  direction  ^  (Figure  3.1),  the  first  term  becomes  a  second  order  derivative  in 
z  alone: 


V2E(z,t)  = 


dz^ 


X  . 


This  leaves  a  differential  equation  in  ^  and  t  which  is  solved  by  separation  of  variables, 


d^E{z,  t)  eoHo  d‘^E{z,  t) 

dz'^  dt"^ 


Because  the  walls  are  perfect  electrical  conductors,  the  transverse  (x)  component  of 
the  electric  field  vanishes  at  2  =  0  and  z  =  L  where  L  is  the  length  of  the  cavity. 
Solutions  are  of  the  form 


k 


jEom  sm{kz)  exp 


rmr 


m 


1,2,3, ••• 


(3.4) 


^Note  here  that  the  hat  denotes  a  vector  direction  and  not  an  abstract  operator  as  exjdained  in 


Section  1.7. 
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with  rn}^  mode  amplitude  ,  propagation  constant  fc,  and  angular  mode  frequency 


VLm  .  Eo^  will  eventually  need  to  be  normalized  when  accounting  for  the  total  energy 
in  the  cavity,  but  for  now  it  will  remain  as  is.  In  order  to  find  the  cavity  mode  fre¬ 
quencies,  this  solution  is  now  substituted  back  into  the  wave  equation.  Equation 
3.3, 

y2p  _  _ 

C2  dt^ 

—  Eq  sin(fc2)  exp^*^’"*^ 

-  {-Ql,Eo  sin(A;2:)  exp’*^"**}  x 

=  Eq  sin(fc2:)  exp”*^”*  x  =  0 

k^  =  (3  5) 


The  frequencies  are  now  defined  according  to 


mCTT 
L  y/e^JIo 


(3.6) 


The  solutions  for  B  are  found  using  E(z,t)  in  Equation  3.5  together  with  Equation 
3.1b, 


V  X  B(z,  t)  = 

dBy{z,  t)  .  _ 

dz  “ 

Bm{z,  t)  = 


epiM)  dE 
c  dt 


.  ^O/^O^tn 

I  - 

C 

.  eo  /.to  m 
*  ck 


■fi'On. 
•  Eq 


m 


sin(A:2)  exp~’^"’‘ 
cos{kz)  exp~*^”' 
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Recall  from  Equation  3.5  that  and  the  expression  for  the  magnetic 

inductance  in  the  cavity  is  obtained, 

Bm{z,t)  =  -  i  y/eo  IM)  cos{kz)  exp“‘^^"'‘  y  .  (3.7) 

3.2  Total  Cavity  Energy 

From  the  last  section,  the  expressions  for  the  mode  field  quantities  £^(-2,  t) 
and  Bm{z,t)  are  shown  below, 

E„,(z,t)  =  FJom  sin(A:„z)  exp“’*^'"*  X 

Bm{z,t)  =  -  i  y/eo  iM)  Eq^  cos{kmz)  exp'’^™*  y  .  (3.8) 

In  general,  the  total  electric  and  magnetic  fields  in  the  cavity  can  be  expanded  in 
terms  of  these  mode  quantities, 


E{z,t) 

=  Eo„  sm{k^z)  exp  x 

(3.9) 

B{z,t) 

=  -  *  V^ofJ-o  cos{kmz)  exp-’^"-*  y  , 

(3.10) 

m 


where  the  amplitude  Eq^  must  now  be  normalized  using  the  total  energy  in  the 
cavity.  The  classical  energy  density  U  in  the  cavity  is  given  by 

lA  =  cq  E  •  E  H - B  •  B 

Using  the  expansions  in  Equation  3.10,  and  averaging  out  the  time  dependence  (the 
(  )  indicates  the  time  averaging),  the  energy  density  becomes 


(ep  /^o) 
2a<o 


■^On 


=  {  f  ^o„,  Eq^,  sin2(A:,„z)  + 

7n,m*  ^ 

■^,  ^  L  ^0 

'ID  'ID  '  ^ 


Integration  of  this  energy  density  over  the  volume  of  the  cavity  produces  the  total 
energy  H  in  the  cavity.  Since  the  field  quantities  are  expanded  in  orthogonal  modes, 
integration  over  the  volume  eliminates  those  quantities  Eq^  Eq^,  for  which  m  ^  m', 

H  =  J  dv i.€o  El^  sm'^{k^z)  +  cos^{k^z)  | 

=  ^  y  du  |eo  sin^(A:,„z)  +  cos^(A:„,z)|  .  (3.11) 


The  volume  of  the  1-D  cavity  is  L  and  the  integral  in  Equation  3.11  is  over  z  to  yield 
the  total  energy, 


(3.12) 


The  expressions  for  E  and  B  in  Equation  3.8  may  be  reexpressed  as 


E,„(z,t)  =  Qmit) 


2nl 


€qL 

BtoC-^)^)  Qmi^^  'v/^OMO 


(m'K  Z\ 

-T) 


X 


^  I'cosI 

frmr 

.epgj  ' 

1  L  ) 

y 


where 


and 


exp 


CqL 


■X  i^m  ^ 


[20?  1 

1 

2 

■  2  ■ 

‘  =  -i  Qrn{t) 

m 

_  _ 

=  Qm{t) 
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The  Qmit)  and  Qm{t)  are  mode  coordinate  and  mode  velocity  variables.  These  new 
expressions  for  E  and  B  are  substituted  into  Equation  3.12  to  yield 


H  = 


Q 


1  L.  . 
+  -2(£o(<o) 


2  +  ‘3m}  • 


(3.13) 


3.3  Translation  to  Quantized  Field 

The  transition  to  a  quantum  mechanical  expression  for  electric  field  as  well 
as  the  total  cavity  energy  involves  a  straightforward  conversion  of  the  and  Ql^ 
variables  to  operators.  Define  Qm  and  Pm  as 

Qm^ 


1 

’  Q=^Ql  =  ^n„Pl 


P  — 

^  m.  — 


Q 


^mQm  —  ^  ^mQ 


^2 


such  that  the  operator  Hp  corresponding  to  the  total  field 


m 


while  the  electric  field  operator  becomes 


E  =  Y^Qm 

m 


sin(fcTO'Z)  • 


(3.15) 


The  raising  and  lowering  operators  for  a  given  field  mode  a^)  corresponding 
to  adding  and  removing,  respectively,  a  photon  from  that  mode  can  be  constructed 
from  the  Q„,  and  Pm  operators: 
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(3.16) 


(^Qm  • 

Solving  these  equations  for  Q.,„  and  Pm  in  terms  of  d„i  and  gives 

Qm  —  (®m  d" 

Pm  —  (^TO  ~  ®m  )  (3.18) 

which  in  turn  can  be  inserted  into  Equation  3.14: 
it,  = 

m 


m 


The  commutator  [d^,  dj„]  is  developed  in  most  quantum  mechanics  texts  and  eval¬ 
uates  to  1,  allowing  the  substitution  d^d^  =  +  1  a>nd 

Hf  =  E^{2aLa™  +  i} 

m 

“  /dj^dro -f- —  1  .  (3.19) 

m  ^  ' 

The  contributed  by  all  modes  leads  to  the  infinite  zero-point  energy  alluded  to 

in  the  first  chapter.  Once  again,  it  is  the  changes  in  energy  with  respect  to  some 

reference  that  can  be  measured,  so  the  problem  can  be  rescaled  to  exclude  the  factor 
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of  I  above.  The  resulting  field  Hamiltonian  is 

Hp  =  Y^ha„,  al„a,n  .  (3.20) 

m 

When  this  operator  acts  upon  a  state  vector  describing  all  modes  of  the  cavity,  the 
total  energy  held  in  each  mode  is  determined  as  the  number  of  photons  n„,  in  the 
mode  times  the  energy  hClm  of  a  single  photon  in  that  mode: 

HF\nin2  •  •  •  Um)  =  {ni  hfli  +n2^f^2  H - H \  n1n2  •  ••  n^)  .  (3.21) 

The  states  labeled  |  ni  n2  •  •  •  n^)  =  |  {^j})  are  in  fact  the  eigenstates  of  the  Hamil¬ 
tonian  Hp- 

The  electric  field  operator  can  now  be  expressed  in  terms  of  the  raising  and 
lowering  operators: 

m 

Meystre  and  Sargent  refer  to  ^  =  £„,  as  the  electric  field  per  photon,  giving 

^  («L  +  «m)  sin(A:„,2:)  .  (3.22) 
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3-4  Matrix  Representation  of  Hp 


The  matrix  elements  of  Hp  are  given  by  |  Hp  |  {ws}),  where 


I  «l)  I  "2)  *  *  ’  l«m> 

<nji  /  nifli  0  •••  0  ^ 

(f>2  I  0  712^2  •  •  •  0 


(nm  I 


0  •  •  '  j 


(3.23) 


By  combining  the  atom  and  field  Hamiltonians,  Ha  and  Hp,  an  uncoupled 
Hamiltonian  can  be  obtained, 

TT  ^  \  '  t 

Ha+F  ~  2  ^  ^  • 

m 

Defining  the  bare  state  representation  as  |  ±  {n^})  =  |  ±)  ®  |  {«,})  yields 


1+00- 

->  |-10...>  • 

1-0-1) 

<+00--.  1 

- 

2 

0 

0 

<-10  -.I 

0 

fil-f  ••• 

0 

a 


"Him  ax 


2 


(3.24) 


This  Hamiltonian  describes  the  energy  of  the  system  when  interactions  between  the 
atom  and  field  are  ignored.  Note  that  since  the  matrix  is  already  diagonal,  the  bare 
states  are  the  eigenstates  of  Ha+p  with  eigenvalues  given  by  the  diagonal  elements. 
For  example,  the  time  independent  Schrodinger  equation  involving  the  |  —  010  •  •  •) 
state  is 

{Ha  +  Hp)\  -olo---)  =  {hn2-^)\  -oio--.). 
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IV.  The  Atom  -  Field  Interaction,  Hint 

4.1  Dipole  Operator  D 

The  classical  energy  of  interaction  between  the  atom  and  field  is  given  by 

Hi„<  =  -er-E  =  D-E,  (4.1) 

where  D  is  the  dipole  moment  of  the  two  level  atom.  The  operator  form  of  the 
equation  above  is  obtained  by  recognizing  the  correspondence  between  the  vector 
quantity  r  and  the  operator  f : 


D  =  — er  D  =  — ef  . 


The  matrix  elements  of  D  define  those  transitions  which  are  dipole  allowed  in  the 
atomic  system.  By  symmetry,  the  matrix  elements  {i\  —  e  r  |  j)  vanish  if  i  and  j 
label  the  same  state.  Since  r  is  antisymmetric, 

(i|-er|j)=  y*  V'*(r)  r  V'(r)  =  0 ,  i  =  j.  (4.2) 

While  it  is  also  possible  for  the  matrix  elements  to  vanish  for  certain  i  7^  j,  this 
model  will  assume  transitions  between  the  two  levels  in  the  atomic  system  are  dipole 
allowed.  The  eigenstates  of  |  ±)  are  used  to  represent  the  dipole  matrix. 


— e 


1+)  I-) 


(+1  /  (+|r|+) 

(+kl- 

<-i  [Hr\+) 

(-kl- 

-e  (  “ 

(+kl- 

1  <-kl+) 

0 
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Let  (+1  —  er  |— )  =  — er+_  =  £>+_  and  (— |  —  er  |+)  =  £>_+  = 


=  L>+_  (7+ +  L>;_  a_  .  (4.3) 

The  matrices  cr+  and  cr_  correspond  to  the  atomic  raising  and  lowering  operators, 
and  d_.  This  can  be  seen  by  their  action  on  the  states  |+)  and  |— )  with  their 
vector  representations,  (10)  and  (01): 


Because  this  system  is  limited  to  two  levels,  raising  of  the  upper  and  lowering  of  the 
lower  states  is  not  allowed.  The  dipole  operator  is  given  by 

D  =  T>+_d+ +  .  (4.4) 


4-2  Interaction  Hamiltonian 

Equations  3.22  and  4.4  together  form  the  quantized  interaction  Hamiltonian, 
Hint-  Combining  the  two  equations. 


D-E 


4-2 


—  (^+  — "1“  ^  ^  (^m  sill{fc/2) 

m 

=  Yli  (^+-^+  +  ^+-^-)  (^m  +  ®m)  sin(A:z)} 

m 

^  XI  ^  {5mat  +  5,„.a^(j+  +  g*r,Mn^-  +  9t,a.m  a^  }  (4.5) 

rn 

where  =  ^-^+-  ’  sin(A;2)  and  g*^  =  sin(fcz)  are  the  coupling 

constant  and  its  complex  conjugate  that  describe  the  strength  with  which  a  given 

atomic  transition  and  field  mode  interact.  With  the  definitions  _ =  — er+_  and 

r  ii/2 

£jn  =  ,  the  explicit  form  of  gm  is 


9m  — 


er+_  cos^ 
2fi 


€oL 


sin(fcz) 


where  cos  9  represents  the  projection  of  the  electric  field  polarization  vector  along 
the  dipole  moment  of  the  atom.  When  the  electric  field  lies  along  the  dipole  axis, 
coupling  is  maximized;  conversely,  a  dipole  aligned  perpendicular  to  the  field  will 
not  couple  to  that  field.  In  the  development  of  the  quantized  field  in  Chapter  3, 
the  polarization  of  the  electric  field  Wcis  taken  to  be  in  the  x  direction  and  as  a 
simplification  the  orientation  of  the  dipole  moment  is  taken  to  lie  in  the  x  direction 
as  well.  The  cos  6  factor  can  then  be  replaced  by  1.  The  coupling  terms  are  then 
given  by 


9m 


er+_ 

2h 


-11/2 


€qL 


sin(fcz) 


(4.6) 


Note  from  Equation  4.5  that  there  axe  four  terms  of  interaction:  9m^m^+, 

g*^alj^a-,  and  g^a.m^-.  The  first  term  describes  a  process  in  which  the  states  of 
the  atom  and  field  mode  are  simultaneously  raised.  Similarly,  the  fourth  term 
describes  the  lowering  of  both  the  atomic  and  field  mode  states.  The  second  and 
third  terms  involve  ‘cross  processes’  in  which  an  exchange  of  energy  occurs  wherein 
the  excitation  of  the  atom  is  accompanied  by  a  lowering  of  the  field  and  vice-versa. 
It  is  shown  in  Meystre  and  Sargent  that  the  contribution  of  the  cross  terms  to  the 
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Hamiltonian  is  on  the  order  of  ^  while  in  comparison  the  contribution  of  the 
remaining  terms  is  [17].  In  general,  the  key  interactions  dealt  with  in  this 

analysis  involve  atomic  and  field  frequencies  which  are  roughly  of  the  same  order  of 
magnitude,  so  the  and  terms  are  neglected.  This  approximation  is 

analogous  to  the  Rotating  Wave  Approximation  in  semi-classical  treatments  of  the 
atom-field  interaction.  The  final  form  of  the  interaction  Hamiltonian  is 

Hint  =  X]  ^  ■  (4-7) 

m 


4-3  Matrix  Representation  of  Hint 

Using  the  properties  of  the  field  and  atom  raising  and  lowering  operators,  the 
bare  state  representation  of  Hint  is  given  by 


1+00 

...)  |-10...> 

1  — j 

(+00-  1 

^  0 

9i 

9m  ^ 

<-10-. •  1 

0 

0 

. 

(4.8) 

<  — -Im  1 

n* 

y  9m 

0 

...  0  J 

The  full  Hamiltonian  Htot  is  given  by 

Htot  “  Ha  +  Hp  +  Hint 

XI  ^  {dm 

9Tn^rn 

1^-}  , 

(4.9) 

m  m 


4-4 


and  the  matrix  elements  of  Htot  are  then  (±  {1^}  |  Htot  |  i  {!«'})  where, 


E 


tot 


(+00...  I 


(-10.  .  I 


(-O.-lm  I 


!+oo...)  1-10..-) 

I 

9*1 


\9 


I  —O'.  -Ini  ) 

9m 

0 

Q  -  - 


2  } 


(4.10) 
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V.  Time  Evolution  of  the  System 


5.1  Time  Evolution  in  the  Bare  State  Basis 

Discussion  of  the  time  evolution  of  the  system  begins  with  the  time-dependent 
Schrodinger  Equation, 

=  (5-1) 

Solution  of  this  equation  for  a  time-independent  Hamiltonian  yields 


I  i>{t))  =  exp 


1^(0)). 


(5.2) 


In  the  bare  state  representation, 


(B  1 7/;{t))  =  {B  I  exp 


i  Htot  t 


\m) 


(5.3) 


where  the  shorthand  notation  \B)  —  |  ±  {la})  has  been  adopted  and  U{t)  = 
exp  |^~*^*‘’**j  is  the  time  evolution  operator.  Using  the  completeness  of  the  bare 
state  basis.  Equation  5.3  can  be  rewritten 


exp 

B' 


(6.4) 


The  matrix  form  of  Htot  is  non-diagonal  in  the  bare  state  basis,  making  the  ex¬ 
ponential  of  the  Hamiltonian  difficult  to  evaluate.  If  the  time  evolution  operator 
is  represented  in  the  eigenbasis  of  Htot,  also  known  as  the  dressed  state  basis,  the 
expression  is  greatly  simplified.  The  completeness  of  the  dressed  states  |  D)  allows 
Equation  5.4  to  be  re-expressed  as 


5-1 


{B\m)  = 


B’  D  D' 


^  H iqi  i 


h 


\D'){D'\B'}{B'\m)  ■ 


(5.5) 


The  exponential  is  expanded  in  a  Taylor  series, 


exp 


Hiftf  t 


h 


h  2n? 


and  the  quantity  {D  \  exp  |  is  evaluated  term  by  term  in  the  expansion. 

The  dressed  states  are  by  definition  orthogonal,  and  the  expansion  terms  evaluate 
as 


(T>|1|T»')  =  1{D\D')=  6d,d' 


(5.6) 


(^1 


%  fitot  ^ 


D') 


i  Ed' t 

h 


{D I  D')  =  —i  lud' t  ^D,D' 


where  Ed'  is  the  eigenvalue  corresponding  to  the  eigenstate  D'  and  ujd'  is  the  associ¬ 
ated  eigenfrequency.  Equation  5.6  can  be  re-summed  to  give  the  resulting  expression 
for  {B\i)(t)), 

(S|V’(t))  =  expf-zwijt]  {D\B'){B'\i){Qi))  .  (5.7) 

B'  D 

There  are  an  infinite  number  of  modes  in  the  cavity  and  therefore  an  infinite  number 
of  bare  states.  The  matrix  form  for  the  time  evolution  in  the  bare  state  basis  is 
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{Bimt)) 


'  {B,\D^)  ■■■ 

e-iu-i  t  Q 

(CilBi) 

(BAm) 

: 

0 

* 

* 

or 

Csit)  =  TUT'^Cb(O)  (5.8) 

where  Cb(<)  are  the  expansion  coefficients  of  ipit)  in  the  bare  state  representation, 


Cs 


'  C+00-{^)  ^ 

^  C_o...l„(t)  y 


(5.9) 


T  is  the  matrix  obtained  by  diagonalizing  Htot  that  transforms  from  the  dressed 
states  representation  to  the  bare  state  representation, 


T 


<+00-| 


|c«l> 

l«»2> 

1  CSm) 

(  Cii 

CI2 

^Im 

C21 

C22 

•  ’  • 

^2m 

^m2 

^mw 

/ 


and  U  is  the  time  evolution  operator  in  the  dressed  state  representation, 

/  exp-‘^i‘/'*  0  •••  ^ 


U 


0  exp 


-i  E2  t/h 


•  •  exp 


0 
0 

-t  Em  t/h 


(5.10) 


(6.11) 
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The  last  two  terms  in  Equation  5.8  are  the  transpose  of  T  and  the  initial  state  of 
the  system  in  the  bare  state  representation. 

The  probability  for  finding  the  system  in  a  given  bare  state  is  found  by  squar¬ 
ing  the  absolute  value  of  the  appropriate  expansion  coefficient.  For  example,  the 
probability  of  finding  the  system  in  a  state  with  an  excited  atom  and  no  photons  in 
the  field  at  time  t  is 


P{+00  •  •  •)  =  <oo...(i)  c+oo-(i)  =  I  c+oo-(i)l''  •  (5.12) 

The  recipe  now  exists  to  compute  the  probabilities  for  the  atom/field  system  to  be 
in  any  given  state  at  a  particular  time.  This  provides  insight  into  how  the  atom 
couples  to  the  various  field  modes  as  a  function  of  time. 


5.2  Truncation  ofMiot 

As  mentioned  in  the  previous  section,  the  infinite  number  of  field  modes  re¬ 
quires  an  infinite  number  of  basis  vectors  to  represent  the  system.  As  a  result,  the 
Hamiltonian  is  infinitely  large.  In  practice,  we  must  decide  how  to  truncate  the 
matrix.  In  other  words,  we  need  to  anticipate  how  many  modes  the  atom  can  be 
expected  to  couple  to.  There  are  two  (related)  ways  to  answer  this.  The  first  is  to 
compare  the  relative  sizes  of  the  coupling  terms  on  the  off  diagonals  to  the  corre¬ 
sponding  term  on  the  diagonal.  For  a  given  value  of  r+_  and  L,  the  coupling  terms 
in  the  first  row  and  column  can  be  written  as  a  constant  times  the  square  root  of  the 
mode  frequency  times  the  spatial  mode  amplitude  modulation  sin(A:z)  =  sin(^^^), 
giving 


where  k  =  The  diagonal  matrix  elements  are 


(jJ 
2  ’ 
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Figure  5.1  Ratio  of  the  off-diagonal  Hamiltonian  terms  to  the  diagonals  as  mode 
number  increases  for  several.  In  this  case  a;  ~  fli. 

and  the  ratio  of  the  off-diagonals  to  the  diagonals  becomes 

sin  [^) 

If  U}  is  comparable  to  the  fundamental  mode  frequency,  and  we  temporarily  ignore 
the  sinusoidal  modulation  in  favor  of  examining  the  envelop  of  the  ratio  values,  then 
the  ratio  varies  approximately  as 

K 

for  large  m.  A  plot  of  this  ratio  is  shown  in  Figure  5.1  for  several  values  of  k. 
If  enough  modes  are  added,  the  off-diagonals  eventually  become  negligible  com¬ 
pared  to  the  diagonals.  As  these  coupling  terms  become  negligible,  the  remainder 
of  the  Hamiltonian  matrix  takes  on  the  appearance  of  the  uncoupled  Hamiltonian 
of  Equation  3.24.  The  Hamiltonian  can  then  be  said  to  be  block-diagonal,  with  a 
coupled  symmetric  block  and  a  diagonal  uncoupled  block  (with  mode  frequencies 
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much  greater  than  the  atomic  frequency).  This  block  diagonal  form  of  is 

^  I  9i  •••  Qm 
9l  •••  0 

[gl  0  ••• 

0 

V 

Hence,  the  atom  is  more  or  less  uncoupled  from  these  higher  modes.  ^  Note  that  if  a 
threshold  is  chosen  below  which  the  coupling  terms  axe  effectively  zero,  then  Figure 
5.1  illustrates  that  larger  couphng  strengths  necessitate  more  modes  and  thus  larger 
matrices.  This  serves  as  the  basis  for  truncating  the  matrix.  Using  the  definitions 
in  Equations  5.9  and  5.10,  the  initial  probability  amplitudes  of  the  eigenstates  (i.e. 
the  expansion  coeflScients)  axe  given  by 


0 

0  ^m+2 

0  0 


a 


0 
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.  (5.13) 
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Ci,(0)  =  T'^CB(0)  (5.14) 
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The  probability  that  the  atom  is  initially  in  the  eigenstate  is  found  by  squaring 
cij.  Recall  from  above  that  by  adding  field  modes,  the  Hamiltonian  can  be  made 
large  enough  that  it  approaches  block-diagonal  form.  The  transformation  matrix 
T  that  diagonalizes  the  Hamiltonian  will  therefore  also  be  block  diagonal.  The 


^It  is  important  to  note  that  regardless  of  the  magnitude  of  uf,  the  matrix  can  always  be  made 
large  enough  to  block  diagonalize  it  by  adding  modes. 
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block  for  which  the  Hamiltonian  is  non-diagonal  will  transform  in  some  complex 
way  giving  rise  to  mixtures  of  eigenstates  comprising  the  bare  states  of  that  block. 
Conversely,  the  diagonal  block  is,  of  course,  diagonalized  by  the  identity  matrix.  The 
transformation  matrix  becomes 
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1  0  •••  0 
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(5.15) 


It  was  shown  above  that  the  distribution  of  eigenstate  probability  amplitudes  for  the 
excited  atom  bare  state  at  t  =  0  was  given  by 


CesM  ^ 

/'on  'I 
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^  ^Iw  j 

A  sample  plot  of  the  initial  eigenstate  probabilities  for  a  coupled  Hamiltonian  in 
which  none  of  the  coupling  terms  were  neglected  is  shown  in  Figure  5.2.  It  confirms 
that  the  eigenstate  amplitudes  do  indeed  eventually  vanish,  so  that  such  plots  as 
these  can  be  used  to  select  an  appropriate  number  of  modes  for  a  particular  coupling 
strength.  In  this  way,  the  amount  of  computation  can  be  reduced  while  obtaining 
the  same  results  (see  Figure  5.3). 
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Probabilities  of  the  eigenstates  comprising  the  initial  excited  atom  state 
of  the  system.  This  system  can  be  truncated  at  roughly  160  modes 
while  still  producing  the  same  result,  as  shown  below. 


Comparison  of  P+oo •••  for  fhe  arbitrarily  large  (200  mode)  and  truncated 
(160  mode)  matrices. 
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5.3  A  Closer  Look  at  Htoi 

Now  that  the  groundwork  has  been  laid  to  perform  the  calculations,  a  closer 
look  at  the  nature  of  the  Hamiltonian,  iI<o<,  should  provide  some  insights  towards 
how  the  system  might  behave  under  several  limiting  circumstances.  The  previous 
section  demonstrated  how  the  Hamiltonian  becomes  block  diagonalized  when  the 
ratio  of  the  couplings  terms  to  the  bare  state  frequencies  (Qm  ~  ^/2)  on  the  diagonal 
vanishes. 

Consider  Figure  5.1  once  again.  As  mentioned  earlier,  the  ratio  vanishes  much 
more  quickly  for  weak  coupling.  It  seems  reasonable  to  assume  that  k  =  — 
be  chosen  such  that  the  atom  couples  to  a  very  small  number  of  modes,  perhaps  even 
one  mode.  In  that  event,  eigenstate  plots  similar  to  Figure  5.2  should  reveal  that 
only  two  eigenstates  (corresponding  to  the  two  dominant  configurations  of  the  sys¬ 
tem)  participate.  For  one  mode  coupling,  the  system  will  exhibit  Jaynes-Cummings 
behavior.  We  shall  see  in  the  next  chapter  that  this  is  in  fact  the  case. 

On  the  other  hand,  when  k  is  significantly  larger,  the  atom  should  couple  to 
many  more  modes.  This  in  turn  implies  that  there  will  be  many  more  eigenstates 
with  non-trivial  amplitudes  arising  from  the  diagonalization  of  •  H  is  intuitive  to 
believe  that  an  atom  coupled  strongly  to  a  large  number  of  modes  should  decay  to  a 
ground  state  with  a  small  likelihood  of  returning  to  the  excited  state.  The  atom  in 
free  space  can  be  considered  the  limiting  case  where  the  atom  couples  to  an  infinite 
number  of  infinitesimally  spaced  modes  of  an  infinitely  long  cavity.  The  second  half 
of  Chapter  6  explores  this  regime  of  strong  coupling  and  narrow  mode  spacing. 

The  next  chapter  details  the  results  of  these  calculations  for  various  values  of 
coupling  strength,  cavity  size,  and  atom-field  detuning.  For  fields  composed  of  many 
modes,  numerical  solutions  are  required.  The  diagonalization  of  the  Hamiltonian  is 
carried  out  by  a  public  domain  linear  algebra  routine  that  returns  eigenvalues  and 
eigenvectors  which  are  then  used  to  propagate  the  system  and  translate  the  result 
to  the  bare  state  basis. 
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VI.  Spontaneous  Emission  Results 

In  this  chapter  the  results  of  the  calculations  based  on  the  previously  developed 
theory  are  presented.  In  the  first  section,  the  one-mode  Jaynes-Cummings  problem 
is  extended  to  the  two-mode  case  to  demonstrate  the  JCM  behavior  of  a  simple 
multimode  system  in  which  the  atom  is  near  resonance  under  weak  coupling  with 
one  of  the  modes.  In  this  regime  of  coupling  strengths,  the  influence  of  the  other 
modes  can  be  neglected.  For  strong  coupling,  the  single-mode  simplification  breaks 
down  and  the  atom  is  shown  to  couple  also  to  the  olF-resonant  mode.  The  field  is 
then  expanded  to  100  modes  under  a  variety  of  coupling  strengths  and  detunings 
to  explore  the  influence  of  the  additional  modes.  Section  6.2  demonstrates  that  an 
atom  accessing  a  large  number  of  closely  spaced  modes  spontaneously  decays  from 
an  excited  state  with  a  rate  comparable  to  an  atom  in  free  space.  All  calculations 
and  results  are  presented  in  atomic  units,  and  conversions  to  MKS  and  Gaussian 
systems  are  provided  in  Appendix  D. 

There  are  essentially  four  parameters  available  to  adjust  in  this  model:  the 
atomic  transition  frequency  a;,  the  position  of  the  two-level  system  z,  the  cavity 
length  L,  and  the  value  assigned  to  In  selecting  values  for  r+_,  it  should  be 
noted  that  it  is  possible  to  choose  a  coupling  strength  that  is  high  enough  to  couple 
the  atom  strongly  to  the  highest  frequency  modes  in  the  field.  One  could  then 
argue  that  there  would  be  influences  from  even  higher  modes  not  accounted  for  in 
the  Hamiltonian.  Such  behavior  is  characterized  by  noticeable  amplitudes  in  the 
probabilities  of  the  outlying  bare  states  as  shown  in  Figure  6.1.  In  order  to  remedy 
this,  either  the  coupling  constant  must  be  reduced  or  the  number  of  modes  must  be 
increased  until  the  outlying  states  no  longer  participate.  In  all  calculations  presented 
here,  the  coupling  was  kept  low  enough  that  a  significant  number  of  the  highest 
states  exhibited  negligible  amplitude  in  the  probability  plots  and  the  amplitudes  of 
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Figure  6.1  Plot  of  all  bare  state  probabilities  for  a  case  in  which  the  coupling 
strength  is  too  high  for  the  size  of  the  Hamiltonian.  The  probability 
of  the  state  corresponding  to  the  excited  atom  is  plotted  at  the  left  on 
the  bare  state  axis.  Note  the  growing  amplitude  of  the  probabilities  for 
the  higher-numbered  bare  states. 
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the  highest  eigenstates  were  negligible  (recall  the  discussion  on  matrix  truncation  in 
Chapter  5). 

One  may  also  adjust  the  coupling  strength  to  each  mode  by  varying  z.  Doing 
so  affects  each  mode  differently  through  the  sm{kz)  term  in  g.  It  follows,  then, 
that  a  particular  choice  of  position  may  cause  one  or  more  of  the  coupling  terms 
to  vanish  if  the  atom  is  located  at  a  node.  For  example.  Figure  6.2  shows  how  the 
choice  of  position  can  affect  coupUng.  In  the  top  diagram,  the  atom  is  located  at 
a  position  where  the  electric  field  amplitude  of  the  fundamental  mode  is  maximum. 
Consequently,  the  coupling  of  the  atom  to  this  mode  is  stronger  than  to  the  second 
cavity  mode  in  (b).  Since  the  amplitude  of  the  second  mode  vanishes  at  the 
coupling  constant  g  =  D+_  ■  £m  sin(A:2)  also  vanishes.  Unless  otherwise  noted,  the 
calculations  presented  below  were  performed  for  a  two  level  system  located  away  from 
the  center  of  the  cavity  at  While  this  has  the  disadvantage  of  modulating  the 
coupling  values  in  a  predictable  yet  uneven  fashion,  it  avoids  the  undesirable  effect  of 
decoupling  half  of  all  of  the  modes  (the  center  position  decouples  the  even-m  modes) 
in  the  cavity. 

6.1  Jaynes-Cummings  Behavior 

As  mentioned  in  Section  1.4,  the  JCM  predicts  an  oscillatory  interchange  of 
probability  between  an  excited  atom  state  and  a  single  excited  field  state  for  which 
the  frequency  of  the  excited  mode  is  closely  tuned  to  the  atomic  transition  (also 
referred  to  as  the  resonant  atom  or  resonant  field  case).  Figure  1.5,  reproduced  here 
as  Figure  6.3  shows  this  complete  interchange.  Now  suppose  that  a  single  mode  is 
added  to  the  existing  field  at  the  next  cavity  mode  frequency  (twice  the  fundamen¬ 
tal  frequency).  When  the  coupling  strength  is  ‘weak’,  an  atom  tuned  near  either 
frequency  will  couple  only  to  that  mode  as  predicted  by  the  JCM.  Figure  6.4  shows 
the  change  in  how  the  atom  is  coupled  to  the  field  modes  as  the  atomic  transition 
frequency  is  ‘dialed’  from  resonance  with  the  fundamental  mode  to  resonance  with 
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Figure  6.2  Strength  of  coupling  for  a  given  mode  is  modulated  by  the  atom’s 
position  in  the  cavity.  Coupling  is  maximized  at  the  antinodes  (a)  and 
zero  at  the  nodes  (b). 


the  second  cavity  mode.  The  atom  is  located  |  of  the  way  along  the  cavity  so  that 
the  sm{kz)  factor  in  the  coupling  constant  contributes  equally  to  both  modes. 

The  Rabi  frequency  can  be  calculated  for  the  two  mode  system  for  comparison 
with  the  analytical  one-mode  result,  where  the  frequency  of  probability  exchange 
went  as  y  =  -f-  4g^.  The  probability  exchange  represented  by  Figure  6.4a  has 

a  half-period  of  T1/2  =  1.865  x  10®  in  atomic  units  (recall  that  the  flopping  frequency 
exhibits  a  cos^  temporal  dependence  in  Equation  1.2,  such  that  a  period  spans  three 
peaks  instead  of  two  for  cos).  This  corresponds  to  an  angular  frequency  of 
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=  3.369  X  10  ®  in  atomic  units  of  frequency. 


For  <5  =  0,  =  2.278  x  lO"®,  (r+_)  =  0.01,  L  =  lOyum  =  1.890  x  10®,  €0  = 

.0795,  e  =  1  and  h  =  1,  all  in  atomic  units,  the  JCM  result  for  y  is 
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Figure  6.3  Characteristic  complete  probability  exchange  predicted  by  the  analyt¬ 
ical  one  mode  JCM  solution  for  a  resonant  atom  (from  Chapter  1). 
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(g)  Atom  is  completely  uncoupled  from  the  two  modes. 
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atom  from  the  first  mode. 


(h)  Equal  detuning  from 
both  modes. 
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Figure  6.4  (a  -  n)  JCM  behavior  of  a  two-mode  system  with  atomic  transition  fre¬ 
quency  tuned  from  resonant  with  fundamental  mode  (a,b)  to  resonant 
with  second  cavity  modes  (m,n).  L  =  lO/itm,  (r+_)  =  0.01. 
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- Two  Mode  Result  o  JCM  Analytical  Result 


Figure  6.5  Overlay  of  the  two  mode  calculation  of  excited  state  probability  with 
the  one  mode  JCM  prediction. 

^  ^  /  2.278  X  10-3  ~ 

2  y  (.0795)  (1.890  X  105) ^ 

=  3.371  X  10“®  in  atomic  units. 


(6.1) 

As  expected,  the  flopping  frequency  predicted  by  the  JCM  is  very  close  to  the  fre¬ 
quency  exhibited  by  the  two  mode  system.  An  overlay  of  the  excited-atom  state 
probabilities  from  the  JCM  and  two  mode  calculations  confirms  the  result  (see  Fig¬ 
ure  6.5).  To  place  the  value  of  the  coupling  strength  in  somewhat  of  a  historical 
context,  the  value  of  er+_  used  in  the  investigation  of  the  ammonia  maser  by  Jaynes 
and  Cummings  [13]  was  1.47  x  10“^^  esu  =  .5787  au  which  corresponds  to  a  coupling 
strength  in  this  scenario  of  gmaser  =  1.95  x  10"'^. 

However,  if  the  coupling  constant  is  significantly  increased,  an  atom  tuned  to 
a  given  mode  can  be  forced  to  couple  to  other  modes  as  well.  This  is  demonstrated 
most  simply  for  the  two  mode  field  by  increasing  the  coupling  by  a  factor  of  200  over 
the  Jaynes-Cummings-like  case  in  Figure  6.4a  Figure  6.6a  indicates  a  significant 

^recall  the  value  of  the  coupling  in  Figure  6.4a  was  r+_  =  .01.  Calculations  indicate  that  r+_ 
values  over  50  are  possible  for  certain  high  quantum-number  states  in  simple  hydrogenic  atoms. 


'Jldata  ^  2  3gg  ^  10“®  au 

=  3.371  X  10"®  au 

2 

2 
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departure  from  the  JCM,  and  complete  breakdown  of  the  one  mode  approximation  is 
achieved  in  6.6b  for  a  factor  of  500  increase  in  coupling.  It  is  clear  that  a  single  mode 
field  is  insufficient  to  adequately  describe  the  spontaneous  emission  properties  of  the 
atom  beyond  the  weakly  coupled  limit.  In  order  to  demonstrate  that  the  presence 
of  many  more  modes  has  a  minimal  impact  on  the  Jaynes-Cummings  behavior  of 
the  system,  we  now  consider  a  field  expanded  to  100  modes  in  which  the  atom  may 
be  slightly  detuned  from  a  mode  near  the  center  of  the  frequency  range,  such  as  the 
49*^  mode.  In  light  of  the  earlier  discussion  of  the  effect  of  atomic  position  on  the 
coupling  terms  in  the  full  Hamiltonian,  the  atom  is  plaeed  at  a  position  selected  to 
decouple  a  minimum  of  modes.  In  this  case,  the  atom  is  placed  ||  x  L  away  from  one 
wall.  Since  the  mode  has  n  + 1  nodes  spaced  ^  x  L  apart,  we  are  assured  that  the 
atom  is  positioned  at  precisely  one  node,  the  19*^  node  of  the  53’’**  mode.  As  before, 
the  cavity  is  small  by  macroscopic  standards,  L  =  10  m/i  =  1.8897  x  10®  a.u.,  and 
the  mode  spacing  is  equivalent  to  the  fundamental  mode  frequency,  2.2278  x  10“^  in 
atomic  frequency  units  (roughly  in  the  infrared).  Once  again,  recall  that  the  effect 
of  detuning  the  atom  from  the  field  mode  in  the  JCM  is  to  reduce  the  amplitude 
of  probability  exchange  and  increase  the  flopping  frequency.  Figure  6.7a  shows  the 
interaction  when  the  atom  is  perfectly  tuned  to  the  mode —  clearly  the  complete 
Hopping  observed  in  the  tuned  JCM.  The  analytical  results  can  be  overlaid  on  the 
solid  line  to  show  how  closely  the  JCM  Hopping  is  duplicated  (see  Figure  6.7b).  Now 
the  atom  is  detuned  from  the  mode  by  5%  of  the  mode  spacing.  The  detuning  causes 
a  reduction  in  the  exchange  amplitude  as  shown  in  Figure  6.8a,  and  the  overlay  of 
the  JCM  prediction  for  this  detuning  in  Figure  6.8b  reveals  excellent  agreement  with 
the  numerical  results.  Now  the  frequency  dial  is  adjusted  to  bring  the  atom  to  10% 
detuning.  The  amplitudes  of  exchange  axe  further  reduced  as  expected,  and  the 
familiar  overlay  of  the  JCM  data  again  indicates  the  agreement  between  analytical 
one  mode  and  numerical  100  mode  results  (Figure  6.9a, b).  As  a  final  check,  the  atom 
is  detuned  to  20%  and  the  probabilities  are  plotted  in  Figures  6.10a  and  6.10b. 
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If  we  now  hold  the  atomic  frequency  fixed  at  10%  detuning  as  in  Figure  6.9, 
we  can  see  how  small  increases  in  coupling  strength  (via  r+_)  can  cause  the  detuned 
atom  to  exchange  more  completely  with  the  field  mode.  In  Figure  6.11,  the  amplitude 
of  the  exchanges  does  indeed  increase  as  r+_  is  raised  from  .05  to  .10,  but  the  overlay 
plot  shows  a  de-phasing  between  the  JCM  prediction  and  the  numerical  results. 
Notice  also  in  Figure  6.12,  where  r+_  =  .15,  that  P+oo  doesn’t  quite  return  to 
one  after  the  initial  decay,  and  the  differences  between  the  analytical  and  numerical 
periods  of  exchange  oscillation  are  increasing.  The  irregularities  that  begin  to  appear 
in  P+oo  a^re  amplified  for  —  .20  (Figure  6.13).  This  is  a  sign  that  the  coupling  is 
strong  enough  for  the  surrounding  modes  to  begin  to  interact  with  the  atom.  Finally, 
a  plot  of  all  of  the  bare  state  probabilities  in  Figure  6.14  confirms  that  this  is  the 
case.  Thus,  in  contrast  to  the  JCM,  where  the  coupling  strength  can  be  increased  ad 
infinitum  until  the  effect  of  detuning  is  negligible  (in  Equations  1.2)  and  exchange  is 
complete,  an  upper  limit  is  placed  on  the  degree  to  which  total  probability  exchange 
can  be  recovered  in  the  multiple  field  mode  cavity. 

The  next  section  investigates  the  other  extreme  of  coupling  strength  and  mode 
spacing,  in  which  a  two  level  system  interacts  with  many  closely  spaced  field  modes 
under  the  influence  of  substantially  higher  coupling  strengths. 
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Figure  6.7  (a,b)  The  atom  is  resonant  with  a  cavity  mode  at  a  frequency  of  .11163 

a.u.  in  a  cavity  with  mode  spacing  .00228  a.u.  In  (b),  the  analytical 
results  from  the  JCM  are  superimposed  on  the  numerical  data. 
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Figure  6.8 


— 100  Mode  Result  o  JCM  Analytical  Result 


(a,b)  Reduction  in  the  amplitudes  of  probability  exchange  result  from 
detuning  the  atom  from  the  field  mode.  Here  the  detuning  is  5%. 
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— 100  Mode  Result  o  JCM  Analytical  Result 


Time  (a.u.) 

(b) 


Figure  6.9  (a)  The  atom  is  detuned  10%  from  the  field  mode.  In  (b),  the  analytical 

results  from  the  JCM  are  superimposed  on  the  numerical  data. 


— 100  Mode  Result  o  JCM  Analytical  Result 


Time  (a.u.) 

(b) 


Figure  6.10  (a)  Detuning  of  20%  of  the  mode  spacing.  Again,  in  (b),  the  analytical 

results  from  the  JCM  are  superimposed  on  the  numerical  data. 
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— 100  Mode  Result  o  JCM  Analytical  Result 


Figure  6.11  (a,b)  The  atom  is  detuned  10%  from  a  field  mode  in  the  cavity,  and  the 

coupling  strength  has  been  doubled  through  an  increase  in  r^_  from  .05 
to  .10.  A  comparison  of  the  JCM  prediction  to  the  numerical  results 
in  (b)  reveals  a  subtle  difference  in  the  frequencies  of  the  probability 
oscillations. 
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Figure  6.12 


— 100  Mode  Result  o  JCM  Analytical  Result 


(a,b)  Further  increase  in  r^._  increases  the  degree  to  which  energy  is 
exchanged  between  the  atom  and  field. 


6-19 


Figure  6.13 
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— 100  Mode  Result  o  JCM  Analytical  Result 


(a)  The  exchange  is  further  increased,  but  the  irregularities  appearing 
in  F+oo  indicate  coupling  to  additional  modes,  (b)  A  non-trivial 
difference  in  the  oscillation  periods  is  quite  apparent. 
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Figure  6.14  The  structure  along  the  first  bare  state  describes  the  evolution  of  the 
excited  atom  state.  Here  the  atom  is  coupling  strongly  to  the  49**^ 
mode  in  the  field,  although  the  ripples  near  the  base  of  the  excited 
mode  probability  indicate  that  coupling  to  a  number  of  off-resonant 
modes  is  non-negligible. 
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Parameter 

Value 

^  Modes 

200 

L 

1.8897x10^  au  =  10  cm 

0,1 

2.2782x10-^  au  =  1.5  GHz 

(Jj 

1.1391x10-5  ay  ^  75  QJJ2 

z 

.6774x10^  au  =  3.58  cm 

Table  6.1  Fixed  parameter  values  for  the  progression  from  weak  coupling  to  strong 
coupling.  r+_  is  varied  to  demonstrate  the  transition. 

6.2  Free  Space- Like  Spontaneous  Emission 

The  results  from  Section  6.1  showed  that  an  atom  tuned  to  a  single  mode  in 
a  multimode  field  exhibited  Jaynes- Cummings-like  behavior  for  systems  with  weak 
couphng  and  widely  spaced  modes.  Here,  the  narrow  mode  spacing  and  strong 
coupling  regime  is  explored.  Figure  6.6  hinted  at  the  fact  that  the  inclusion  of 
additional  modes  for  stronger  coupling  can  produce  markedly  different  behavior  than 
is  predicted  by  the  JCM.  In  fact,  under  the  influence  of  strong  coupling  with  a  large 
number  of  modes,  the  atom  exhibits  spontaneous  emission  rates  similar  to  those  of 
an  excited  atom  in  free  space.  The  excited  atom  state  probability  is  shown  to  follow 
a  decay  which  is  approximately  exponential  in  the  high  coupling  regime,  and  the 
rate  of  this  decay  is  compared  to  the  rate  of  decay  predicted  for  an  atom  in  a  one 
dimensional  cavity  by  the  Wigner-Weisskopf  theory  [17]. 

To  demonstrate  this,  we  begin  with  an  uncoupled  system  in  which  the  atom  is 
resonant  with  a  mode  of  the  cavity — in  this  case,  the  50*^  cavity  mode — but  g  =  0. 
Imagine  now  that  there  is  a  dial  controlling  the  value  of  with  which  the  atom-field 
interaction  may  be  adjusted.  This  dial  has  four  settings:  the  initial  zero  coupling 
r+_  =  0.0  and  three  higher  settings  of  r+_  =  0.225,  r+_  =  0.450,  and  r+_  =  4.500. 
By  progressing  through  these  settings,  we  can  see  how  the  behavior  of  the  system 
evolves  towards  a  nearly  free  space  decay  of  excited  atom  state  probability.  The 
fixed  parameters  (number  of  modes,  L,  f2i,  w,  and  z)  for  this  series  of  scenarios  are 
summarized  in  Table  6.1.  In  Figure  6.15a,  the  atom  is  completely  uncoupled  from 
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the  field.  As  a  result,  the  atom  remains  in  the  excited  state.  The  accompanying 
plot  of  the  eigenstate  amplitudes  shows  that  the  system  is  described  entirely  by 
one  eigenstate,  consistent  with  the  uncoupled  scenario.  In  the  next  two  figures,  the 
coupling  has  been  turned  up  so  that  the  atom  interacts  not  only  with  the  resonant 
mode,  but  also  with  several  off-resonant  modes  as  well.  The  off-resonant  coupling  is 
indicated  by  the  increasingly  rapid  and  erratic  probability  evolutions.  Finally,  the 
last  figure  in  the  series.  Figure  6.18,  demonstrates  what  appears  to  be  an  exponential 
decay  in  the  |  -t-  000  •  •  •)  probability.  When  all  bare  state  probabilities  are  plotted 
as  in  Figure  6.19,  the  extent  of  participation  by  the  large  number  of  modes  becomes 
evident. 

The  decay  from  this  multimode  analysis  is  compared  in  Figure  6.20  to  the 
exponential  decay  predicted  by  the  Wigner-Weisskopf  Theory.  The  presentation 
of  this  theory  in  Meystre  and  Saigent  is  based  on  a  three  dimensional  vacuum  of 
sufficient  size  that  the  modes  form  a  continuum.  When  the  theory  is  appropriately 
reduced  to  the  one  dimensional  case  through  the  application  of  the  1-D  mode  density 
and  the  limitations  on  the  orientations  of  the  dipole  moment  and  field  polarization, 
the  predicted  decay  rate  F  is 

er+_ 
ceok 

such  that  the  excited  atom  state  probability  falls  as  exp~^*. 

Using  the  values  of  u)  and  r+_  that  produced  Figure  6.18a,  the  one  dimen¬ 
sional  Wigner-Weisskopf  theory  predicts  an  exponential  free  space  decay  rate  of 
r  =  2.115x10“®  in  inverse  atomic  units  of  time  (8.744  x  10^^  s“^),  while  an  expo¬ 
nential  fit  to  the  results  in  Figure  6.18  gave  a  decay  rate  of  3.70  x  10“®  in  inverse 
atomic  units  of  time  (1.53  x  10^^  s“^).  Figure  6.18a  is  reproduced  in  Figure  6.20  with 
an  overlay  of  the  Wigner-Weisskopf  prediction  for  comparison.  However,  it  must  be 
stressed  that  the  behavior  of  the  two  level  system  in  the  closely  spaced  mode  cavity 
is  not  irreversible  as  in  the  free  space  spontaneous  decay.  The  theory  of  Wigner  and 
Weisskopf  predicts  a  strictly  exponential  decay  and  does  not  allow  the  possibility  of 
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revivals  in  the  excited  atom  state  probability  after  a  finite  period  of  time.  In  con¬ 
trast,  these  revivals  do  occur  for  the  quantized  system  presented  here.  Figure  6.21 
illustrates  this  probability  revival  that  occurs  when  the  temporal  axis  of  Figure  6.18a 
is  extended  by  a  factor  of  100.  The  revivals  are  a  result  of  the  phasing  properties  of 
the  eigenstates.  Since  each  eigenstate  evolves  at  its  own  eigenfrequency,  then  it  is 
reasonable  to  expect  that  there  exists  times  at  which  groups  of  the  eigenstates  are 
in  phase  (recall  that  at  t  =  0  all  eigenstates  are  in  phase),  causing  resurgences  in  the 
amplitude  of  the  initial  bare  state  of  the  system. 
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Probability  (g  ProbabiUty 


(a)  Excited  Atom  State  Probability 


(b)  Eigenstate  Amplitudes 


6.15  I  +  000  •  •  •)  probability  and  eigenstate  amplitudes  for  r+_  =  0.0.  The 
atom  is  completely  uncoupled  from  the  field  and  remains  in  the  excited 
state. 


(a)  Excited  Atom  State  Probability  (b)  Eigenstate  Amplitudes 


Figure  6.16  |  +  OOO***)  probability  and  eigenstate  amplitudes  for  r+_  =  0.225. 

Coupling  to  several  off-resonant  modes  introduces  irregularities  in  the 
time  evolution. 
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(a)  Excited  Atom  State  Probability  (b)  Eigenstate  Amplitudes 


Figure  6.17  |  +  000  •  •  •)  probability  and  eigenstate  amplitudes  for  r+_  =  0.45. 


(a)  Excited  Atom  State  Probability  (b)  Eigenstate  Amplitudes 


Figure  6.18  |  +  000  •  •  •)  probability  and  eigenstate  amplitudes  for  r+_  =  4.50. 


Bare  State 


Figure  6.19  Evolution  of  the  bare  state  probabilities  for  r+_  =  4.5.  The  exponen¬ 
tial  decay  of  the  excited  atom  state  results  from  the  strong  coupling 
with  many  modes. 


_ Multimode  Analysis  a  Wigner-Weisskopf  Result 


Figure  6.20  Comparison  of  Wigner-Weisskopf  spontaneous  emission  decay  and  re¬ 
sults  from  Figure  6.18a. 


Figure  6.21  Revivals  in  the  excited  atom  state  probability  appear  when  the  system 
is  allowed  to  evolve  for  much  longer  timescales  (note  the  change  in  the 
Time  axis  from  Figure  6.18). 
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—  L  =  10  cm 


- - L  =  5  cm 


Figure  6.22  The  atom  decays  from  its  excited  state  more  slowly  in  the  smaller 
cavity  in  which  the  mode  spacing  is  larger. 

Now  consider  the  effect  of  reducing  the  cavity  size  on  the  decay  rate.  As  the 
cavity  size  is  reduced,  the  mode  spacing  becomes  larger.  An  atom  tuned  near  one 
mode  of  the  cavity  becomes  more  detuned  from  the  rest  of  the  modes.  Therefore,  if 
u  and  r+_  are  held  fixed  while  L  is  reduced,  then  the  atom  will  couple  less  to  the 
off-resonant  modes  and  the  decay  rate  should  decrease.  Figure  6.22  compares  the 
decays  of  a  cavity  of  length  L  =  10cm  =  1.8897  x  10^  a.u.  (solid  line)  and  the  cavity 
reduced  by  half  to  L  -  5cm  =  9.4486  x  10®  a.u.  (dashed  line).  The  two  level  system 
was  still  positioned  at  ||  L.  Even  for  such  a  relatively  small  change  in  cavity  length, 
the  decay  is  noticeably  more  gradual  for  the  smaller  cavity. 

6.3  Unanticipated  Effects 

The  results  of  the  previous  two  sections  illustrated  that  appropriate  choices 
of  the  various  parameters  lead  to  the  recovery  of  results  for  two  limiting  cases  in 
agreement  with  two  independent  well-established  theories  of  spontaneous  emission. 
However,  the  appropriate  choices  did  not  reveal  themselves  easily.  In  exploring 
the  various  combinations  of  coupling  strengths,  atomic  transition  frequencies,  cavity 
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lengths,  and  matrix  truncation  points,  several  regimes  were  discovered  in  which  the 
behavior  of  the  coupled  atom- field  system  could  be  characterized  as  anomalous.  In 
the  first  case,  a  two-level  system  exhibiting  complete  Rabi  flopping  with  a  mode 
could  be  caused  to  decouple  from  the  mode  by  drastically  increasing  the  number  of 
modes  in  the  field.  As  the  mode  number  increased,  the  amplitude  of  the  probability 
exchanges  decreased  as  if  the  atom  were  being  detuned  from  the  field.  The  same 
effect  also  manifested  itself  in  the  investigation  of  free  space-like  behavior  discussed 
in  Section  6.2:  the  exponential  spontaneous  decay  illustrated  in  Figure  6.18  became 
more  gradual  and  eventually  evolved  to  a  probability  curve  which  dropped  quickly 
from  1  to  a  steady  state  value  near  .7  when  the  number  of  modes  was  expanded  to 
1000.  It  seems  likely  that  this  may  be  a  result  of  an  ill-conditioned  Hamiltonian 
matrix,  in  spite  of  the  fact  that  the  diagonalization  routine  indicated  convergence 
of  the  solutions  and  the  conservation  of  probability  and  energy  were  verified  (see 
Appendix  A).  The  first  step  in  addressing  this  problem  will  be  to  investigate  the 
stability  of  the  Hamiltonian  under  a  variety  of  coupling  strengths  and  matrix  sizes, 
and  analyze  in  depth  the  method  by  which  the  diagonalization  is  carried  out  by  the 
LAPACK  routine.  If  numerical  instability  in  the  matrix  is  ruled  out  as  a  cause  for 
this  behavior,  then  the  conclusion  which  remciins  is  that  this  is  simply  the  way  the 
system  behaves  under  those  conditions,  and  the  ‘problem’  so  to  speak  lies  in  the 
interpretation  of  these  results. 

A  similar  phenomenon  occurred  for  the  Jaynes-Cummings-like  system  in  which 
the  atom  was  perfectly  tuned  near  a  single  cavity  mode.  When  the  perfectly  tuned 
100  mode  system  in  Section  6.1  was  subjected  to  gradually  increasing  coupling 
strengths  through  increases  in  r+_,  the  spontaneous  emission  was  again  suppressed. 
An  increase  of  a  factor  of  ten  in  r^_  resulted  in  a  5%  reduction  in  the  exchange 
amplitude.  However,  recall  from  Section  6.1  (Figures  6.11-6.13)  that  the  detuned 
atom  behaved  in  the  opposite  fashion — the  exchanges  in  the  probability  plots  became 
more  complete  as  r^._  was  increased.  Several  plots  illustrating  the  decoupling  effect 
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on  the  perfectly  tuned  atom  are  included  in  Appendix  C,  but  no  further  explanation 
is  provided  here  for  this  behavior. 
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VII.  Conclusion 


The  scale  of  today’s  microlaser  devices  are  quickly  approaching  sub-micron  dimen¬ 
sions  as  military  and  commercial  requirements  push  the  envelop  of  fabrication  tech¬ 
nologies.  Microcavity  structures  such  as  the  VCSEL  oflFer  the  promise  of  large  arrays 
of  highly  efficient  lasers  to  meet  a  variety  of  needs.  The  forward  progress  in  making 
these  microcavity  lasers  more  efficient  inevitably  relies  on  a  thorough  understanding 
of  the  spontaneous  emission  processes  within.  This  thesis  develops  a  fully  quantum 
mechanical  model  of  the  spontaneous  interaction  between  a  two-level  atom  com¬ 
prising  the  gain  medium  and  a  multimode  electromagnetic  field  within  an  optical 
resonator  cavity.  The  research  described  in  the  previous  chapters  provides  a  foun¬ 
dation  for  future  exploration  into  the  behavior  of 

7.1  Quantization  of  the  Atom,  Field,  and  Their  Interaction 

One  of  the  first  steps  in  addressing  a  problem  quantum  mechanically  is  to 
develop  a  Hamiltonian  for  the  system  under  observation.  In  this  case  the  system 
consists  of  a  two-level  atom  interacting  with  a  multimode  electromagnetic  field  within 
a  one  dimensional  perfectly  conducting  cavity.  As  a  simplification,  only  processes 
in  which  a  single  photon  was  either  emitted  from  or  absorbed  by  the  atom  were 
considered.  Chapter  2  describes  the  steps  taken  to  construct  the  atomic  Hamiltonian. 
The  result  is  given  in  Equation  2.7  as 


where  is  the  Pauli  z  spin  operator.  A  representation  in  the  orthogonal  two  level 
atomic  basis  {  |  +) ,  |  — )  }  can  be  taken  in  which  the  matrix  elements  of  Ha  are 

+  \Ha\  +  )  = 
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In  Chapter  3,  the  classical  wave  equation  is  solved  subject  to  the  boundary  conditions 
for  a  PEC  cavity.  The  expression  for  the  electric  and  magnetic  fields  are  normalized 
to  the  total  electromagnetic  energy  in  the  cavity  and  the  field  is  quantized  by  asso¬ 
ciating  each  mode  with  a  quantum  mechanically  treated  harmonic  oscillator.  The 
Hamiltonian  Hp  for  the  field  is 

Hp  —  h  Cljn  ^In^m 

m 

with  photon  raising  and  lowering  operators  and  (Equation  3.20).  The  basis 
with  which  to  represent  Hp  consists  of  states  in  which  there  are  no  photons  in  any 
of  the  modes  or  one  photon  in  one  of  the  modes:  {  |  0000  •••),  |0---lm0---)}.  The 
matrix  elements  of  Hp  are 

(OOOO--- |/fF|0000---)  =  0 
(liOO---  |i7ip|li00---)  = 

(00---l™|ifir|00..-U>  =  hflm 

Since  Ha  and  Hp  commute,  the  uncoupled  operator  Ha+f  can  be  constructed: 

hu)  ^  V  >  t  ^ 

Ha+f  2  ^  ^ 

m 

with  an  eigenbasis  formed  from  the  direct  product  of  the  atomic  and  field  bases 
{  I  +  {0})  5  I  ~  {I3})  }•  These  eigenstates  of  the  uncoupled  Hamiltonian  Ha+f  are 
the  bare  states  of  the  system  and  are  useful  in  describing  the  state  of  the  system  be¬ 
fore  and  after  transitions  have  taken  place.  Chapter  4  treats  the  interaction  between 
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the  atom  and  the  field.  The  interaction  is  given  by 

Hint  —  ^  “I" 

7ri 

where  d+  and  d_  are  the  atomic  raising  and  lowering  operators  and  p„i  and  its 
conjugate  characterize  the  strength  of  coupling  between  the  atom  and  mode 
of  the  field.  In  the  bare  state  representation,  Hint  couples  the  excited  atom  state, 
I  +  {0}))  fo  each  of  the  states  |  —  {1^}  ).  The  full  Hamiltonian  is  therefore  given  in 
Equation  4.9  as 

Htot  =  Ha  +  Hp  +  Hint 

TYi  m 

The  matrix  representation  of  Htot  in  the  bare  state  basis  is  non-diagonal  when  the 
coupling  terms  are  non-zero,  so  it  is  useful  to  find  a  basis  in  which  the  Hamiltonian 
is  diagonal.  The  transformation  matrix  T  that  diagonalizes  the  matrix  IHIfoi  consists 
of  elements  which  axe  the  projections  of  the  eigenbases,  or  dressed  states,  onto  the 
bare  states.  The  time  evolution  of  the  bare  state  expansion  coeflScients  of  the  system 
wavefunction  is  then  computed  using  Equation  5.8: 

Csit)  =  TUT‘Cs(0) 

The  probability  to  find  the  system  in  a  particular  bare  state  at  time  t  is  obtained 
by  squaring  the  absolute  value  of  its  expansion  coelRcient  in  •  The  behavior  of 
the  system  is  analyzed  by  inspecting  the  way  in  which  these  probabilities  evolve  in 
time.  The  results  of  these  analyses  obtained  in  Chapter  6  for  the  limiting  cases  of  an 
atom  weakly  coupled  to  a  single  field  mode  and  an  atom  strongly  coupled  to  many 
modes  are  summarized  below.  Following  Section  7.2,  recommendations  are  provided 
for  future  studies  beyond  the  scope  of  this  project. 
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1.2  Overview  of  the  Results 

7.2.1  Jaynes- Cummings  Limit.  In  the  limit  of  weak  coupling,  an  atom 
in  a  two  mode  field  was  shown  to  behave  in  a  manner  consistent  with  the  Jaynes- 
Cummings  Model,  an  analytical  model  which  describes  the  interaction  of  a  two  level 
atomic  system  with  a  single  field  mode.  The  presence  of  additional  modes  was  shown 
to  have  minimal  impact  for  an  atom  near  resonance  with  and  weakly  coupled  to  one 
mode,  but  increases  in  the  coupling  strengths  through  r+_  caused  the  atom  to  begin 
coupling  to  off-resonant  modes. 

7.2.2  Free  Space- Like  Behavior.  When  the  two  level  atom  was  placed  in 
a  field  of  many  tightly  spaced  modes  and  the  dipole  matrix  element  proportional 
to  was  increased  by  a  factor  of  ~100  over  those  used  in  the  Jaynes-Cummings 
limit,  the  atom  coupled  to  a  large  number  of  modes  and  decayed  exponentially  to 
the  lower  state.  The  decay  rate  exhibited  by  the  atom  was  comparable  to  the  rate 
calculated  using  the  Wigner-Weisskopf  Theory  of  spontaneous  emission  for  an  atom 
in  free  space;  3.70  x  10“®  in  inverse  atomic  units  of  time  (1.53  x  10^^  s“^)  from  a  fit 
of  the  data  compared  to  the  Wigner-Weisskopf  result  of  2.115x10“®  in  inverse  atomic 
units  of  time  (8.744  x  10^^  s“^).  Investigation  of  the  effect  of  changing  the  cavity  size 
indicated  that  reduction  in  L  caused  a  decrease  in  decay  rate,  a  result  attributable 
to  the  increase  in  mode  spacing  and  therefore  decrease  in  the  effective  coupling  of 
the  atom  the  modes  in  the  cavity. 

7.2.3  Anomalous  Results.  In  the  pursuit  of  these  investigations,  several 
regimes  of  parameter  values  produced  results  which  were  counterintuitive  and  with¬ 
out  explanation.  In  one  regime,  it  was  discovered  that  the  addition  of  a  significant 
number  (500  -  800)  of  modes  to  a  system  which  clearly  exhibited  either  Jaynes- 
Cummings  or  free  space-like  behavior  caused  that  behavior  to  change  drastically. 
In  both  cases,  the  atom  in  the  presence  of  the  expanded  field  tended  to  remain  al¬ 
most  to  the  point  of  exclusion  in  the  excited  state.  It  can  be  seen  by  looking  at  the 
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form  of  the  Hamiltonian  matrix  that  it  is  possible  to  make  the  matrix  large  enough 
through  the  addition  of  modes,  and  hence,  states,  that  the  ratio  of  the  largest  to 
smallest  elements  can  become  quite  large  (~  10^  —  lO'^).  This  indicates  a  rather 
ill-conditioned  matrix  which  may  have  caused  numerical  instabilities  in  the  diago- 
nalization  of  the  Hamiltonian.  Without  an  in-depth  analysis  of  the  diagonalization 
methods  used  by  the  LAPACK  routine  employed,  it  is  difficult  to  draw  any  concrete 
conclusions  regarding  possible  connections  between  the  anomalous  results  and  the 
condition  number  of  the  system.  In  another  regime,  a  moderate  increase  in  cou¬ 
pling  strengths  could  cause  a  suppression  of  decays  out  of  the  excited  atomic  state 
for  an  atom  perfectly  tuned  to  a  cavity  mode  and  previously  exhibiting  complete 
probability  exchanges.  As  shown  in  Section  6.1,  however,  the  detuned  atom  did 
not  experience  this  suppression.  This  phenomenon  is  not  easily  attributed  to  the 
conditioning  of  the  matrix,  and  no  explanation  can  be  offered  here. 

7.3  Recommendations  for  Future  Study 

As  alluded  to  in  the  abstract,  the  research  detailed  in  this  thesis  carries  an 
exploratory  element  to  it.  Along  the  way  I’ve  encountered  more  than  a  few  dead 
ends,  numerous  dead  ends  that  turned  out  not-so-dead,  and  a  great  many  questions 
as  yet  without  answers.  It  is  in  that  context  that  I  provide  my  recommendations  for 
possible  future  studies  for  which  this  research  may  in  part  serve  as  a  foundation. 

1.  The  anomalies  referred  to  in  the  previous  section  need  to  be  investigated  in 
depth.  While  it  is  quite  possible  that  numerical  instabilities  are  in  part  to 
blame  for  the  behavior  of  the  large  matrices,  the  inhibition  of  spontaneous 
decay  due  to  coupling  increases  for  a  perfectly  tuned  atom  manifests  itself  in 
what  appear  to  be  stable  systems. 

2.  In  this  analysis,  the  simplification  was  made  to  treat  the  cavity  as  a  perfect 
conductor.  However,  the  cavities  within  VCSELs  and  other  microcavity  lasers 
are  by  no  means  perfect  conductors.  A  natural  complement  to  this  research  is 
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an  exploration  of  the  imperfectly  conducting,  or  lossy,  cavity.  In  more  realis¬ 
tic  treatments  of  microcavities,  there  is  a  certain  penetration  of  the  radiation 
into  the  walls  of  the  resonator,  and  for  all  practical  purposes  it  is  necessary 
to  artificially  introduce  loss  at  one  end  to  obtain  an  output  for  the  photons 
oscillating  within.  One  key  implication  of  a  lossy  cavity  is  the  introduction 
of  a  finite  linewidth  to  the  emitted  radiation.  The  details  of  integrating  this 
finewidth  into  the  Hamiltonian  are  to  be  worked  out,  but  it  may  be  plausible 
to  assume  that  each  cavity  mode  with  its  inherent  lineshape  can  be  discretized 
into  a  collection  of  tightly  packed  modes  of  infinitesimal  hnewidths  whose  am¬ 
plitudes  follow  a  Lorentzian  shape.  The  behavior  of  an  atom  coupled  to  a 
single  longitudinal  mode  could  be  ascertained  by  examining  the  interaction  of 
the  atom  with  the  modes  under  the  Lorentzian  profile. 

3.  A  related  oflFshoot  addresses  the  problem  of  the  dispersive  cavity.  Instead  of 
the  atom  in  a  vacuum  cavity,  the  gain  medium  should  be  modeled  to  more 
accurately  reflect  the  characteristic  of  the  semiconductor  materials  used  in 
laser  diodes.  This  has  several  implications.  First,  the  behavior  of  the  medium 
itself  must  be  addressed.  Secondly,  the  dispersive  nature  of  the  semiconductor 
materials  which  fill  the  cavity  has  an  impact  on  the  mode  structure  of  the 
electromagnetic  field.  The  two  level  system  modeled  here  more  accurately 
describes  the  quantum  dot  structures  which  are  yet  another  area  of  current 
interest  in  microcavity  edge  emitting  and  surface  emitting  lasers. 

4.  In  its  infancy,  this  research  was  intended  to  explore  the  three  dimensional  cubic 
cavity.  In  order  to  simplify  some  of  the  issues  related  to  the  coupling  terms, 
the  problem  evolved  into  the  one  dimensional  cavity  addressed  here. 

All  of  these  recommendations  are  based  on  a  need  to  actively  pursue  an  under¬ 
standing  of  the  spontaneous  processes  in  microcavity  devices.  This  thesis  constitutes 
a  first  step  towards  addressing  the  interactions  between  the  gain  medium  and  the 
multimode  field  that  influence  laser  efficiencies  at  an  atomic  level.  Ultimately,  this 
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research  complements  the  existing  body  of  knowledge  employed  by  Air  Force  sci¬ 
entists  and  Research  Planners  in  improving  current  microcavity  technologies  and 
assessing  the  direction  of  tomorrow’s  laser  programs. 
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Appendix  A.  Verification  of  the  FORTRAN  90  Code 

One  of  the  most  important  steps  in  carrying  out  a  numerical  analysis  of  a  complex 
problem  such  as  this  is  to  verify  that  the  computer  code  used  produces  valid  results. 
This  is,  of  course,  a  misleading  statement — the  analyses  are  usually  employed  to 
address  problems  for  which  ‘valid’  results  are  the  unknowns.  However,  it  is  typi¬ 
cally  the  case  that  exploratory  research  is  built  upon  prior  results  in  which  there  is 
sufficiently  high  confidence.  This  appendix  describes  the  steps  that  were  taken  to 
validate  the  FORTRAN  90  code  in  Appendix  B  to  as  high  a  degree  as  possible.  The 
process  of  validation  involved  certain  checks  described  briefly  as  follows: 

Does  the  code  compute  the  cavity  mode  frequencies  correctly? 

The  mode  frequencies  were  computed  manually  on  a  PC  using  a  simple  16-bit  pre¬ 
cision  function  created  in  Mathematica  for  a  variety  of  cavity  lengths.  These  were 
used  as  a  baseline  for  comparison  with  the  results  from  the  code.  The  purpose  was 
to  catch  gross  errors  in  the  algorithms  versus  compare  the  numerical  accuracy  of  the 
machines  used.  Quite  obviously,  the  capabilities  of  the  SGI  Origin  2000  for  produc¬ 
ing  high  accuracy,  high  precision  results  heavily  outweigh  the  Pentium  66.  However, 
a  comparison  showed  that  the  results  agreed  to  seven  or  eight  significant  digits. 

Does  the  code  construct  the  full  Hamiltonian  correctly? 

Once  again,  the  objective  was  to  catch  errors  in  the  way  the  matrix  elements  were 
computed.  Sample  matrices  were  generated  by  another  Mathematica  routine  and 
compared  to  the  results  from  the  code,  then  both  of  these  were  checked  against 
a  matrix  generated  by  hand.  In  matrices  as  large  as  100x100,  the  diagonals  and 
coupling  terms  sampled  at  random  compared  as  favorably  with  the  code  as  the 
frequencies  above. 
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Are  the  results  from  the  diagonalization  routine  reliable? 

To  verify  the  performance  of  the  LAPACK  diagonalization  routine  and  its  interface 
with  my  code,  the  eigenvalues  and  eigenvectors  of  1,  2,  5,  20,  and  100  mode  sys¬ 
tems  were  compared  with  those  generated  by  the  built-in  Eigensystem  function  in 
Mathematica.  An  additional  self-consistency  check  is  provided  in  the  code;  for  each 
run,  a  routine  computes  the  product  of  the  eigenvector  matrix  with  its  transpose. 
Since  these  matrices  perform  unitary  transformations  of  the  Hamiltonian,  the  prod¬ 
uct  T'^T  should  result  in  the  identity  matrix.  One  of  the  outputs  of  the  code  is 
a  sum  of  the  elements  of  this  product  which  in  all  cases  should  produce  the  total 
number  of  columns  or  rows  {an  n  xn  identity  matrix  logically  has  n  entries  of  1  for 
its  diagonals).  This  sum  was  confirmed  following  each  run. 

Does  the  code  evolve  the  system  in  time  properly? 

Part  of  this  verification  is  provided  as  the  body  of  Section  6.1.  The  code  results  for 
1,  2,  5,  and  20  mode  systems  were  overlaid  on  the  analytical  (for  the  1  mode  system) 
and  Mathematica  results  with  no  discernible  differences  in  the  plots. 

Is  probability  conserved  for  the  system  as  it  evolves  in  time? 

It  is  a  basic  necessity  that  at  all  points  in  time  the  total  probability  for  observing  the 
system  in  one  of  its  states  is  1.  In  order  to  confirm  this,  the  probability  amplitudes 
for  all  states  must  be  squared  and  summed  for  all  times  t.  One  of  the  routines  in 
the  code  calculates  these  sums  for  every  time  sampled  and  writes  them  to  a  file.  In 
all  cases,  the  sum  of  the  probabilities  was  1  to  within  5  significant  digits. 

Is  energy  conserved  for  the  system  as  it  evolves  in  time? 

Since  there  is  no  loss  mechanism  built  into  the  system,  it  is  expected  that  a  cavity 
which  contains  a  specified  total  energy  at  time  t  =  0  should  have  that  same  total 
energy  at  each  time  step  thereafter.  A  routine  was  built  into  the  code  to  compute 
the  total  energy  in  the  system  for  each  iteration  over  the  time  variable.  The  results 
of  this  were  then  written  to  yet  another  file.  Because  the  analyses  all  started  with 
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the  system  in  the  excited  atom  state,  the  initial  energy  in  the  cavity  was  the  energy 
of  the  excited  atomic  state:  The  results  written  to  the  energy  conservation 

file  consistently  demonstrated  at  every  time  step  that  the  total  energy  in  the  cavity 
remained 

Can  the  original  full  Hamiltonian  be  reconstructed  using  the  eigenvectors  and 
eigenvalues  from  the  diagonalizing  routine? 

The  most  compelling  test  of  the  diagonalizing  routine  is  the  reconstruction  of  the 
original  full  Hamiltonian  from  the  diagonal  Hamiltonian  and  the  transformation 
matrix  both  obtained  from  the  LAPACK  routine.  A  comparison  of  the  reconstructed 
and  original  Hamiltonians  provides  a  figure  of  merit  T  in  assessing  the  accuracy  of 
the  LAPACK  solutions.  The  reconstructed  Hamiltonian  is  found  from 

itoi  =  TTE' 

where  IF  is  the  diagonalized  dressed  state  (diagonal)  representation  of  the  Hamilto¬ 
nian.  IHIfot  can  then  be  subtracted  element-wise  from  the  differences  squared, 
then  the  square  root  of  the  sum  of  the  resulting  squares  computed,  and  finally  nor¬ 
malized  to  the  matrix  size  by  dividing  by  n^.  For  the  a  sampling  of  the  matrices 
evaluated  in  this  thesis,  the  values  of  T  fell  around  In  other  words,  the  orig¬ 

inal  and  reconstructed  Hamiltonians  agreed  element  by  element  out  to  around  the 
tenth  or  eleventh  digit. 
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Appendix  B.  Source  Code 


Main  Program: 

This  program  computes  the  time  evolution  of  the  bare  states  for  a 
quantized  coupled  atom-field  system  in  a  one -dimensional  perfect  electrically 
conducting  cavity.  The  atomic  system  is  quantized  as  a  2-level  system  and  the 
field  is  comprised  of  from  1  to  n  normal  modes,  where  n  can  be  as  large  as 
computational  resources  will  allow.  A  coupled  Hamiltonian  is  created  based 
on  user  input  parameters,  then  is  diagonalized  using  a  public  domain 
LAPACK  routine.  The  eigenvalues  and  eigenvectors  axe  then  used  with 
the  user-specified  initial  state  of  the  system  to  determine  the  time 
evolution  in  terms  of  the  bare  states.  The  components  are: 

decs.module  -  Module  containing  the  global  declarations 
funct .module  -  Module  containing  several  common  functions 
readinl.ld  -  Subroutine  to  road  the  user  specified  parameter  values 
fullham.ld  -  Subroutine  to  construct  the  full  Hamiltonian  and 

diagonalize  using  the  LAPACK  routine  dsyevd.for 
timodev.ld  -  Subroutine  to  compute  the  time  evolution  of  the 
system  in  the  bcLre  state  basis 


PROGRAM  mainlD 
USE  decs. module 

USE  funct.raodule 
IMPLICIT  NONE 


CALL  readin.ld 
CALL  fullham.ld 
CALL  timedov.ld 

END 


B-l 


Modules : 

The  following  two  modules  are  structures  which  contain  information 
or  functions  which  are  global  entities.  The  module  decs.module 
contains  the  global  declarations,  while  funct_module  contains 
two  functions  used  by  several  subroutines. 


MODULE  decs_module 
IMPLICIT  NONE 
Constants: 


DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 


PRECISION  c 
PRECISION  pi 
PRECISION  epsO 
PRECISION  muO 
PRECISION  me 
PRECISION  hbar 
PRECISION  e 


Values  from  Bransden&Joachain 


PARAMETER  (c  =  137.O36D0  ) 
PARAMETER  (pi  =  3 . 14169265359D0  ) 
PARAMETER  (hbar  =  l.ODO  ) 
PARAMETER  (epsO  -  0. 07967747 155D0  ) 
PARAMETER  (muO  =  6 .69176247674D-4  ) 
PARAMETER  (e  =  l.ODO  ) 
PARAMETER  (me  =  l.ODO  ) 


Array-related  declarations : 


DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 
DOUBLE  PRECISION,  ALLOCATABLE  :: 


INTEGER 
INTEGER 
INTEGER 
CHARACTER (2) 
CHARACTER(2) 
CHARACTER (2) 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 


imax 

iterate 

count max 

amethod 

fmethod 

gmethod 

gstart 

testval 

Iworkfactor 


a_f_array(: , :) 
iutarray( : , : ) 
fullarray( : , : ) 
tempf ull ( ; , ; ) 
afreqarray( : ) 
f f reqarray ( : ) 
phot_num€Lrray( : ) 
gvector( :) 
eigenvalues ( : ) 
bareinitiaK : ) 


1 /alpha 


l/4Pi 

4Pi/c-‘2 


Atomic  system  dec Icurat ions ; 

DOUBLE  PRECISION  rba 
DOUBLE  PRECISION  a.freqhalf 
DOUBLE  PRECISION  z 
Field  declarations: 


INTEGER  m 
INTEGER  mmax 
INTEGER  phot^num 
DOUBLE  PRECISION  1 

Initial  conditions : 

INTEGER  i state  ! 

Time -development  declarations : 

DOUBLE  PRECISION  tmax 
DOUBLE  PRECISION  tmin 
DOUBLE  PRECISION  tstep 

END  MODULE  decs.module 


mode  number 
total  number  of  modes 
number  of  photons 
cavity  dimensions 


specifies  the  initial  state  of  atom 


Module  containing  common  functions  for  field  frequencies 
and  mode  spatial  distributions 

MODULE  funct .module 

USE  decs.module 
IMPLICIT  NONE 


CONTAINS 

Fiinction  to  calculate  the  radiation  frequencies 

DOUBLE  PRECISION  FUNCTION  r.freq  (mode, dim) 

INTEGER  mode  !  mode  numbers 
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DOUBLE  PRECISION  dim 
r_freq  =  c*pi*mode/dim 


cavity  dimensions 


END  FUNCTION 

Function  to  calculate  the  value  of  the  field  spatial 
distribution  at  a  given  location  for  a  given  mode 

DOUBLE  PRECISION  FUNCTION  g^spatial (mode, dim, pos) 

INTEGER  mode  !  mode  numb 

DOUBLE  PRECISION  dim  !  cavity  dii 

DOUBLE  PRECISION  pos  !  atom  posi 


mode  numbers 

cavity  dimensions 

atom  position  in  cavity 


g_spatial  =  SIN( 
RETURN 

END  FUNCTION 
END  MODULE  funct_module 


SIN(mode  *  pi  ♦  pos  /  dim) 
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Subroutine  readin^ld: 

Opens  a  stream  to  the  file  ’microcav .dat*  and  reads 
the  user  input,  assigning  values  to  the  appropriate  variables. 
The  conversions  from  MKS  to  atomic  units  occur  in  this 
subroutine  where  necessary. 


SUBROUTINE  readin^ld 


USE  decs_module 
USE  funct^module 
IMPLICIT  NONE 

INTEGER  i ,  j 
CHARACTER(l)  fgheader 
CHARACTER(13)  dummy 
CHARACTER (8)  day 
CHARACTER (10)  clock 
DOUBLE  PRECISION  Imicrons 
in  microns 

DOUBLE  PRECISION  zmicrons 
DOUBLE  PRECISION  atomic 
DOUBLE  PRECISION  fund.freq 
DOUBLE  PRECISION  step 
DOUBLE  PRECISION  terapwave 
DOUBLE  PRECISION  tempg 


!  declarations 
!  common  functions 

!  iterators 

!  header  letter  at  top  of  f^or.g.inp  file 
!  dummy  for  reading  in  (skipping  over)  strings 
!  used  for  date  stamp  in  message  file 
!  used  for  time  stamp  in  message  file 
!  accepts  input  cavity  length 

!  input  position  in  microns 

!  input  transition  wavelength/frequency  depending  on  amethod 
!  sets  fundamental  mode  frequency  if  f method  is  manual 
!  used  for  creating  the  gvector  array  of  coupling  constants 
!  temporary  holder  for  read-in  field  wave lengths /frequencies 
!  temporcLry  holder  for  read-in  g  values 


Open  a  message  file  for  run-time  error  statements  and  other 
miscellaneous  junk  (time/date  stamped) 


OPEN  (unit=7,  f i le= * me s sages. dat * ,status=* replace*) 
CALL  DATE_ AND.TIME (D ATE=day , TIHE=clock) 

WRITE  (7,  1)  day. clock 
WRITE  (7,  >(‘'Messagos:”)*) 


Open  the  input  file 


OPEN  (unites,  file=*micro_ld. inp* ,status=*old* ) 
Read  in  the  values 


READ (8, 15)  amethod 
READ (8,15)  f method 
READ (8. 15)  gmethod 
READ (8, 30)  gstart 
READ (8. 10)  dummy 
READ (8, 30)  atomic 
READ (8, 30)  rba 
READ (8. 10)  dummy 
READ (8. 30)  Imicrons 
READ (8. 30)  zmicrons 
READ (8. 20)  mmar 
RE AD (8. 20)  phot_num 
READ (8, 10)  dummy 
READ (8. 20)  iterate 
READ(8.10)  dummy 
READ (8. 20)  istate 
READ (8. 35)  tmin 
RE AD (8. 35)  tmax 
READ (8. 35)  tstep 
READ (8, 35)  testval 
READ (8. 20)  Iworkfactor 


Initialization  of  input -dependent  quantities 

imax  =  mmax+1  !  array  dimensions 

Convert  from  input  microns  to  the  a.u.  used  internally 

l=lmicrons*18897 . 26877774 
z=zmi crons* 18897 . 26877774 


Flag  to  make  sure  2-level  system  is  in  the  cavity! 

IF(z.GT.l  .AND.  gmethod.EQ. *a*)  THEN 

WRITE (7,  *  ("micro. Id. in  Error  —  Verify  z<l’')*) 


Array  initialization  block 

ALLOCATE ( intaxray ( imax , imax) ) 
ALLOCATE ( a_  f .array ( imax , imax ) ) 
ALLOCATE(fullarray (imax, imax) ) 
ALLOC  ATE  (tempfulK  imax,  imax)  ) 
ALLOCATE (afreqarray( iterate) ) 
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AL L 0 C ATE ( f f re qarr ay ( mmax ) ) 

ALLOCATE (eigenvalues (imax) ) 

ALLOCATE (bareinitial(imax) ) 

ALLOCATE(phot_nmnarray(iinax) ) 

ALLOCATE (gvectorC imax) ) 
intarray  =0.0 

a^f^array  =0.0  !  matrices 

fullarray  =  0.0 

eigenvalues  =0.0  !  vectors 

bareinitial  =  0.0 

This  block  calculates  atomic  frequencies  based  on  input  file  values 
SELECT  CASECeunethod) 

CASE( *  AW* , *awO  !  Single  frequency  is  calculated  based  on  atomic  wavelength 
!  The  factor  0.0455633  converts  from  micron  wavelengths 
!  to  a.u.  frequencies 
af reqarray (1 )  =  0 . 0455633 /atomic 

CASE( * AF* , *af * )  i  A  single  frequency  is  read  in  directly  from  micro_ld.inp 

af reqarray (1)  =  atomic 


CASE( *M* , *m* )  !  Creates  array  of  length  *  iterate*  of  atomic  frequencies 

!  The  values  range  from  0  to  twice  the  fundamental  frequency 
fund_freq  =  r_freq<l,l) 
af  reqeurray  (1 )  =f  und^f  req 
IF(iterate.GT.l)  THEN 

af reqarray ( 1 ) =0 . 0 

step  =  2*fund.freq/(iterate-l) 

DO  i=2, iterate 

af reqarray ( i ) =af reqarray ( i - 1 ) +st ep 

END  DO 

END  IF 

END  SELECT 


Create  an  array  of  field  frequencies  using  either  the  r.freq  function 
or  read-in  wavelength  values  from  f_or_g.in 

ff  reqarray  =  O.OdO 
SELECT  CASE(fmethod) 

CASE( *AA* , *aa* )  !  calculate  field  resonant  frequencies 

DO  i=l,mmax 

ffreqarrayd)  =  r_freq(i,l) 

END  DO 

CASE(*HW* , *mw*)  !  read  in  field  wavelengths 


IF(gmethod.EQ.*al*  .OR.  gmethod.EQ .  *m* )  THEN 

OPEN (unit=10 , action=  *  READ  * , f ile=  *  f _or_g . inp  * , status=  *  old  * ) 
READ (10, 15)  fgheader 
IF(fgheader.EQ. *f *)  THEN 


ELSE 


tempwave=l .  OdO 
DO  i=l,mmax 

READ (10, 30)  tempwave 
ffreqarrayd)  ®  0 . 0465633/terrmwave 

END  DO 

WRITE(7 ,  *  ('*f_or_g .  in  —  Verify  header  matches  data”)*) 
STOP 


END  IF 


ELSE 


WRITE(7 ,  *  ("micro^ld.in  —  gmethod  must  not  be  & 

&  mf  {i.e.  read  wavelengths  OR  g  values  from  file,  not  both}”)*) 
STOP 


END  IF 

CASE(*MF* , *mf *) 


read  in  field  frequencies 


IF(gmethod.EQ.  *al*  .OR.  gmethod.EQ.  *m*)  THEN 

OPEN (unit=10 , action=  *  READ  * , f ile=  *  f _or_g . inp  * , status=  *  old  * ) 
READ (10, 15)  fgheader 
IF (fgheader , EQ . *  f  * )  THEN 
tempwave=l . OdO 
DO  i=l,mmax 

READ (10, 30)  ffreqarrayd) 

END  DO 


ELSE 


WRITE(7 , d”f_or  V. in  —  Verify  header  matches  data”)*) 
STOP 


END  IF 


ELSE 


WRITE(7,  *  (”micro„ld.in  —  gmethod  must  not  be  & 
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END  IF 


&  mf  {i.e.  read  wavelengths  OR  g  values  from  file,  not  both}")*) 
STOP 


1 

5 

10 

15 

20 

30 

35 

40 


END  SELECT 


The  €Lrray  phot.numarray  holds  the  number  of  photons  in  the  field 
phot^numarray  =  phot^num 

This  block  initializes  the  array  holding  the  values  of  g  for  each 
mode  which  are  artificially  reduced  for  increasing  mode  number 
arithmetically  from  gstart  down  to  0.  Note  that  gvector(l)  is  always 
set  to  0  since  this  corresponds  to  a  dipole  forbidden  +/+  atomic 
transition 


SELECT  CASE(gmethod) 

CASE(*M1 * , *ml* , *  A* , *a* , *A2* , *a2*)  !  makes  all  g  values  =  1.0;  if  *a* ,  gvector 

!  will  be  overwritten  in  sub  fullham  Id 

gvector  =  1 , OdO 
gvector(l)  =  O.OdO 

CASE( ’MF* , *mf * )  !  read  in  g  values 

gvector (1)  =  0.0 

IF (f method. EQ . *AA*  .OR.  fmethod.EQ . *aa* )  THEN 

OPEN  (unit=10 ,  actions^  *  READ  * ,  f  ile=  *  f  _or_g .  inp  * ,  status=  *  old  * ) 

READ (10, 15)  fgheader 
IF (fgheader . EQ . *  g  * )  THEN 
DO  i=l,nimax 

READ (10, 30)  tempg 
gvoctor(i+l)  =  tempg 

END  DO 

ELSE 

WRITE(7, * ("f_or_g,in  —  Verify  header  matches  data")*) 
STOP 

END  IF 

WRITE(7,  *  ("micro_ld.in  --  gmethod  must  not  be  & 

&  mf  *(i.e.  read  wavelengths  OR  g  values  from  file,  not  both}")*) 
STOP 


ELSE 

END  IF 

CASE(*A1* ,*al*) 


!  calculate  g  values  from  field  frequencies 

gvector (1)  =0.0 
DO  i=l,  mmax 

gvector (i+1)  =  SQRT(ffreqarray(i) ) 

END  DO 


END  SELECT 


Initialize  the  bare  state  expansion  coefficients 
IF(istate.GT.imax)  THEN 

WRITE (7,  *  ("micro ..Id.  in  Error  —  Verify  istate  <-fmmax+2}")  * ) 
STOP 


ELSE 

bareinitiaK  istate)  =  1.0 

END  IF 


F0RMAT(* Date/Time:  *,a9,lX,alO) 

FORMAT(all) 

FORMAT(a) 

FORMAT (a5) 

F0RMAT(i6) 

F0RMAT(fl6.1O) 

FORMAT (f 20. 10) 

F0RMAT(a8) 

RETURN 

end  !  Subroutine  readin.ld 
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Subroutine  fullham.ld: 

This  subroutine  populates  the  full  Hamiltonian  then  diagonalizes  it  using  the 
LAPACK  routine  DSYEVD.  The  eigenvectors  and  eigenvalues  are  saved  to 
’combined.dat’  for  later  use.  ^  ^  _ _ _ 


SUBROUTINE  fullham.ld 

USE  decs^module  !  declarations 

USE  funct^module  !  common  functions 

IMPLICIT  NONE 

iterators 

temporary  mode  numbers 
constants  in  g 

spatial /mode-dependent  component  of  g 
field  frequency 
temporarily  holds  eigenvalues 


Define  several  functions  for  use  in  populating  the  arrays: 

element (m) :  combines  the  field  frequency  computation 

with  the  spatial  dependence  of  the  matrix 
element 

rad^freqCm) :  calculates  field  frequency  using  spec’d  dimensions 

element (m)  =  SQRT(r_freq(m,l) ) 
rad_freq(m)  =  r_freq(m,l) 

Initialize  LAPACK  parameters 

Iwork  =  (l+6*imax+2*imax*(INT<L0G10(REAL(imax))/LQG10(2.0))  +  1)  +  3*imax**2)  *  Iworkf actor 
liwork  =  2  +  6*imax 

ALLOCATE(work(lwork) ) 

ALLOCATE (iwork (liwork) ) 

work  =0.0 
iwork  =  0 


INTEGER  i ,  j ,  count 
INTEGER  mtemp 
DOUBLE  PRECISION  g.const 
DOUBLE  PRECISION  element 
DOUBLE  PRECISION  rad.freq 

DOUBLE  PRECISION,  ALLOCATABLE  ::  evarray(:,:) 

LAPACK  declarations 

INTEGER  Iwork 
INTEGER  liwork 
INTEGER  info 

INTEGER,  ALLOCATABLE  ::  iwork(:) 

DOUBLE  PRECISION.  ALLOCATABLE  ::  work(:) 


INTERACTION  HAMILTONIAN 

This  section  populates  the  interaction  hamiltonian  using  the 
coupling  term  g  =  -  (e  rba  /  2  hbar)  [4  hbar  Wn  /  V  e0]''l/2  f(x,y,z) 
if  using  gmethod  =  *a*  (automatic)  or  forces  g*s  to  specified  values 
if  gmethod  =  'm’ . 

Open  a  series  of  files  for  output  data  including  uncoupled  hamiltonian 
terms,  interaction  terms,  dressed  state  eigenvalues,  and  the  full 
hamiltonian 

OPEN  (unit=9 ,  f  ile=  *  combined .  dat  * ,  status=  *  replace  * ) 

OPEN (unit=12 , f ile=  *  f ullham . dat ’ , st atus=  *  replace  * ) 

Determine  whether  to  calculate  or  force  g  values  and  carry  it  out 
SELECT  CASE  (gmethod) 

CASE(»Ml^»ml^»MF^»mf»,»Al^*alO  !  Sets  Hint(l,i)  =  g(i) 

DO  i=l ,  imax 

intarray ( 1 , i )  =  gvector(i)  *  SQRT(DBLE(phot_num) ) 

END  DO 

CASE(*A\'aO  !  Calculates  g(i)  forHint(l,i) 

g_const  =  -  e  ♦  rba  ♦  SQRT(  hbar/(l*epsO)  )  /  (2*hbar) 
intarrayd  ,1)  =  0.0 
gvector(l)  =  0.0 
m  =  1 

DO  i=2,  imax  !  this  step  overwrites  initial  gvector  values 

mtemp  =  m 

gvector (i)  =  g_const  *  element (mtemp)  *  g.spatial (mtemp, l,z) 
intarrayd  ,i)  =  gvector (i)  ♦  SQRT(DBLE(phot_num)  ) 
m  =  mtemp+1 

END  DO 

CASE(  *  A2  *  ,  *a2*  )  !  Calculates  g(i)  for  Hintd.i) 

g_const  =  -  e  *  rba  ♦  sqRT(  hbar/ (l*epsO)  )  /  (2*hbar) 
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intarray (1 ,1)  =0.0 
gvector=  0.0 

gvector(2)  =  g.const  *  element (1)  *  g_spatial(l , 1,2) 

DO  i=2,  imax 

intarray (1  ,i)  =  gvector(i)  *  SQRT(DBLE(phot_num) ) 

END  DO 

END  SELECT 

First  column  is  same  as  first  row 
DO  i=2,  imax 

intarray (i ,1)  =  intarray (1 ,i) 

END  DO 

Interaction  hamiltonian  created 


Repeat  the  calculation  countmax  times  to  process  each  atomic  frequency 
specified 

This  block  determines  how  many  times  (countmax)  to  repeat  the  calculation 


SELECT  CASE(amethod) 

CASE(»AF*  ,»af  *  ,»AW’,  »awO 
countmax=l 
CASE(*M*,*mO 

countmax=iterate 


END  SELECT 


ALLOCATE  (  evar  ray  (  imax ,  countmax)  ) 

DO  count=l , countmax 

t  DO  BLOCK 


ATOM-FIELD  UNCOUPLED  HAMILTONIAN 

The  following  block  populates  the  atom-field  uncoupled  array 

a_f_array(l  ,1)  =  afreqarray(  count) /2  !  +  (phot_num-l)*rad_freq(mtemp) 

m  =  1 

DO  i=2,  imax 

mtemp  =  m 

a_f .array (i,i)  =  -afreqarray( count) /2  +  phot _num*f freqar ray (i-1) 
m  =  mtemp+1 

END  DO 

Uncoupled  hamiltonian  created 


FULL  HAMILTONIAN: 


Simple  statement  to  add  the  uncoupled  and  interaction  Hamiltonians; 

Since  the  input  array  to  dseyvd  is  altered  to  hold  the  eigenvectors, 
create  a  temporary  full  hamiltonian  to  pass  to  the  diagonalizing  routine 


fullarray  =  hbar  *  (a_f .array  +  inteirray) 
DO  i=l,imax 

If (fullarray ( 1 , i) /fullarray (i , i ) 
fullarray ( 1 , i ) =0 . 0 
fullarray ( i , 1 ) =0 . 0 

END  IF 


END  DO 


.LT.  testval)  THEN 


tempfull  =  fullarray 
Diagonalize 

CALL  dsyevd(’V* , *U' 1 imax, tempfull, imax, eigenvalues, work, 1 work, iwork,liwork, info) 


Write  the  resulting  eigenvalues  to  a  file  and  append  for  each  iteration  over 
the  atomic  frequencies 

WRITE(7,  * ("Diagonalizing  routine  output:  ")*) 

WRITE (7,  1)  info 

WRITE(9,5)  count,  afreqar ray (count) 

DO  i=l,imax 

WRITE(9,20)  (tempfulKi,  j)  ,  j=l,imax) 

END  DO 

WRITE (9, 20)  (eigenvalues (i ) ,  i=l,imax) 

DO  i=l,imax 

e varray ( i , count ) =e igenvalues ( i ) 

END  DO 
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Write  the  full  and  interaction  hamiltonians  to  a  file  and  append 
for  each  iteration 


\^i^miiititi*:tf.**^***3it*********m**************  end  do  block  ★*********«*s»i*****»**+****3f«*>tt**:4i***)** 
END  DO 

!  Write  results  to  appropriate  files 

!  DO  i=l,iinax 

!  WRITE (11, 40)  (evarray(i, j) ,  j=l , countmax) 

!  END  DO 

DO  i=l , imax 

WRITE(12,50)  (fullarrayd, j)  ,  j  =  l,imax) 

END  DO 


1  F0RMAT(i4) 

5  F0RMAT(i4,2X,el2.6) 

10  FORMATC’  100000(e9.3,  IX)) 

20  FORMATdOOOOO  (el8 . 10 ,  IX)  ) 

30  F0RMAT(*  *,  100000(el2 .6,  IX)) 

40  F0RMAT(»  100000(el3 .7,  IX)) 

50  F0RMAT(»  \  100000(el2 .6,  IX)) 

! print  *,  ^pre-deallocation  okay* 

DEALLOCATE(iwork) 

DEALLOCATE(work) 

DEALLOCATE (a_f_array) 

DEALLOCATE ( intarray) 

!  DEALLOCATE(fullarray) 

D£ALLOCATE(afreqarray) 

DEALLOCATE(ffreqarray) 

DEALL  OC  ATE  (phot  _nuinarray  ) 

DEALLOCATE  (gvector) 

DEALLOCATE (eigenvalues ) 

! print  *,  *fullhain  completed* 

RETURN 

END  !  Subroutine  fullham_ld 
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Subroutine  timedev.ld: 

Uses  the  eigenvectors  and  eigenvalues  generated  in  fullhamld  to 
evolve  the  system  in  time  and  traoislate  the  result  to  the 
bare  state  basis  _  _  _ 


SUBROUTINE  timedev.ld 

USE  decs .module 
IMPLICIT  NONE 


INTEGER  i ,  j  ,  count , 

INTEGER  dummyl 
DOUBLE  PRECISION  duinmy2 
DOUBLE  PRECISION  t 
DOUBLE  PRECISION  eigencheck 


timeiteration 


C0MPLEX(8),  ALLOCATABLE 
C0MPLEX(8),  ALLOCATABLE 
C0MPLEX(8),  ALLOCATABLE 
C0MPLEX(8),  ALLOCATABLE 
DOUBLE  PRECISION,  ALLOCATABLE 
DOUBLE  PRECISION,  ALLOCATABLE 
DOUBLE  PRECISION,  ALLOCATABLE 
DOUBLE  PRECISION,  ALLOCATABLE 


Uii(:,:) 

:  eigen ( : , : ) 

:  eigentranspose ( : , : ) 
:  time_evolve( :) 

:  probeirray  (  : ) 

:  Gvprodarray( ; 
:  hamtest( : , 

:  newhamC : , : ) 


) 


Allocate  space  for  the  arrays 

ALLOCATE (Uii (imax ,  imax) ) 

ALLOCATE(eigen(imax,imax) ) 

ALLOCATE (eigentranspose (imax , imax) ) 

ALLOCATE (time_evolve( imax) ) 

ALLOCATE (probarray( imax) ) 

ALLOCATE (evprodarray ( imax ,  imax ) ) 

ALL OCATE(nevham( imax, imax) ) 

ALLOCATE (hamtest ( imax , imax) ) 

OPEN (unit=13 , f ile=  *  timedev .dat ’ , status=  *  replace  * ) 

OPEN (uni t=14 , f ile=  *f irstate . dat  * , status=  *  replace  * ) 

OPEN (\mit=15,file=’ probsq.dat* ,status= 'replace  * ) 

0PEN(\init=16,file='conserv .dat '  ,status=' replace') 

OPEN (unit=17 , f ile= ' evproduct . dat ' , status= ' replace ' ) 

OPEN (unit=18 , f ile= 'hamtest . dat ' , status= 'replace ' ) 

OPEN (unit=19 , file= 'realtime . dat  * , status= ' replace ' ) 

OPEN (\mit=20 , f ile= ' imtime . dat ' , status= ' replace ' ) 

REWIND (9)  !  resets  the  read  pointer  back  to  the  beginning  of  the  file 

Initialize  some  arrays 

Uii  =  CMPLX(0. 0,0.0) 
time.evolve  =  CMPLX(0 .0,0 .0) 

Iterate  the  time -evolution  computation  for  each  of  the  atomic  frequencies 
DO  count=l ,  countmax 


BLOCK******************’*'*********’*'*********** 


Read  the  eigenvalues  and  vectors  from  the  file  'combined.dat' 

READ (9, 5)  dummy l,dummy2 

'tempfull'  now  becomes  a  temporary  holder  for  eigenvectors  then 
eigenvalues 

Read  eigenvectors 

tempfull  =  O.OdO 
DO  i=l,  imax 

READ (9, 20)  (tempfull (i ,j) ,  j=l,imax) 

END  DO 

evprodarray  =  MATMUL (tempfull, TRANSPOSE ( terapf ull ) ) 

DO  i=l,  imax 

WRITE(17,20)  (evprodarray(i , j) ,  j=l,imax) 

END  DO 

The  following  matrices  are  the  eigenvector  matrix  and  its  transpose 
for  calculating  the  transition  probabilities 

eigen  =  CMPLX( tempfull) 
eigentranspose=  TRANSPOSE (eigen) 

Verify  that  the  eigenvectors  are  orthomormal 

eigencheck  =  REAL (SUM (MATMUL (tempf ull, TRANSPOSE(tempf ull)))) 

WRITE(7, '  (**Sum  of  the  elements  of  eigenvectors  x  (eigenvectors)T : ") ' ) 
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WRITE (7, 30)  eigencheck 
!  Read  eigenvalues 

tempfull  =  O.OdO 

READ (9, 20)  (tempfull(l,i) ,  i=l,imax) 

!  Write  the  eigenvalues  back  to  *  comb ined.dat*  to  verify  they ^ re  read  correctly 

WRITE(9,20)  (tempfull(l,i) ,  i=l,imax) 

!  Verify  the  diagonalization  is  accurate  by  multiplying  Ev .  Uii  .  EvT  to  recover  Hf ull 

newham=0 .0 
DO  i=l,imax 

newhajn(i)i)  =  tempfull  (1  ,i) 

END  DO 

hamtest  =  MATMUL ( REAL ( eigen)  ,  MATMUL(newham,REAL(eigentranspose) )) 

DO  i=l,imax 

WRITE(18,20)  (hamtest(i, j) ,  j=l,imax) 

END  DO 

!  Now  the  fun. . .for  t=0  to  tmax  in  steps  of  tstep,  the  vector  of  time- 

!  dependent  bare-state  expansion  coefficients  is  calculated  from 

i  I  I  Ci(t)  I  I  =  I  I  Eigenvec  I  I  .  M  Uii  11.11  EigenvecT  | {  .  | |  Ci(0)  | | 

timeiteration  =  1 
t  =  tmin 

DO  WHILE(t.LE.tmax) 

DO  i=l,  imax 

Uii(i,i)  =  CMPLX(COS(terapfull(l,i)  *  t) , -SIN (tempfull (l,i)  *  t)  ,8) 

END  DO 

time^evolve=MATMUL(eigen,MATMUL(  Uii,MATMUL(eig6ntranspose,bareinitial) )) 

!  Write  real  and  imaginary  parts  of  time.evolve  to  files  for  diagnostic  purposes 

WRITE<19,26)  t,  REAL(time.evolve(l)) 

WRITE(20,25)  t,  AIMAG(time_evolve(l)) 

!  Compute  probabilities 

DO  i=l,imax 

prob€u:ray(i)  =  REAL(CONJG(time_evolvo(i) )  ♦  time^evolve(i) ,8) 

END  DO 

!  Write  the  results  for  the  current  value  of  t  to  files: 

WRITE(13,25)  t,  (probarray(i) ,  i=l,imax) 

WRITE (14, 25)  t,  probarray(l) 

!  Verify  conservation  of  probability  and  energy 

WRITE(16,25)  t,  SUM(probarray) 

WRITE ( 16 , 25)  t ,  DOT.PRODUCT  (t ime_evolve , MATMUL  (f ullarray ,  t ime_evolve)  ) 

timeiteration  =  timeiteration  +  1 
t  =  t  +  tstep 

END  DO  !  end  of  time  iterations 

END  DO  !  end  of  atomic  frequency  iterations 


Notes: 

Uii  contains  the  terms 


1 

Exp[-i  El  t/hbar] 

0  0 

1 

0  Exp[-i 

E2  t/hbar]  0 

Uii(i,i)  =  1 

0 

0  ExpC-i  E3  t/hbar] 

1 

1 

1 

0 

where  El,  E2,  are  the  eigenvalues 


!  Format  statements 

5  F0RMAT(i4,2X,el2.6) 

10  FORMAT(a) 

20  FORMAT(100000(el8 . 10 , IX) ) 

25  F0RMAT(el4 . 6 , lx , 100000 (el8 . 10 , IX) ) 

30  F0RMAT(fl0.5) 

RETURN 

END  !  Subroutine  timedev_ld 
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Appendix  C.  Inhibited  Spontaneous  Emission  Plots 

For  the  sake  of  completeness,  several  plots  demonstrating  the  suppression  of  decay 
out  of  the  excited  atomic  state  alluded  to  in  Chapter  6  are  provided  on  the  following 
pages.  No  additional  discussion  of  the  phenomenon  is  provided  at  this  time. 
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Probability 


Time  (a.u.) 

Figure  C.2  Same  system  for  r+_  =  0.05. 
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Probability  Probability 


-  P |+00->  “  “  ^1-10-) 


-  P |+00->  -  -  P l-io-.) 


Figure  C.4  =  0.5. 
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-  P I+00-)  -  -  P I-10  -) 
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Appendix  D.  Unit  Conversions 


Listed  below  are  several  conversions  for  converting  between  atomic  units  and  other 
unit  systems.  Values  are  taken  from  Bransden  and  Joachain  Appendix  11  [7]  as  well 
as  the  Naval  Research  Laboratory  Plasma  Formulary  [6]. 


Quantity 

Atomic  Units 

MKS 

Gaussian 

Angular  Momentum 

h 

1.0546  X  10"^^  J  s 

1.0546  X  10-21  erg  s 

Charge 

e 

1.6022  X  10-^®  C 

4.8032  X  Ifl-i®  stat.coul 

Length 

«o 

5.2918  X  10"“  m 

5.2918  X  10-®  cm 

Velocity 

^0 

2.1877  X  10®  m  s-^ 

2.1877  X  10*  cm  s-i 

Time 

Ofl. 

Vo 

2.4189  X  10-^^  s 

2.4189  X  10-11  g 

Frequency 

2  TT  ao 

6.5797  X  10^®  s-i 

6.5797  X  lOi®  s-i 

Energy 

4  TT  €o  OO 

4.3598  X  10-1®  j 

4.3598  X  10-11 

Note  that  the  value  of  en  in  atomic  units  is  A-. 
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Following  is  a  graphical  representation  of  a  portion  of  the  electromagnetic 
spectrum  to  provide  a  scale  for  comparison  of  atomic  iinit  frequencies  to  more  familiar 
quantities. 


Region 

A 


Microwave 


X  nm 


V  Hz 


50,000  5.0x10 


12 


13 


to  au 
.00076 


- 

-  10,000 

3.0x10 

.00456 

Infrared 

14 

1  ^  800 

3.8x10 

.05775 

Visible 

= 

14 

400 

7.5x10 

.11399 

Ultraviolet 

Figure  D.l  A  portion  of  the  electromagnetic  spectrum  for  comparison  of  atomic 
units  of  frequency  to  transition  wavelengths  (nm)  and  transition  fre¬ 
quencies  (Hz). 


D-2 


Bibliography 


1.  Abdel-Hafez,  A.  M.,  A.-S.  F.  Obada  and  M.  M.  A.  Ahmad.  “N-Level  Atom  and 
N-1  Modes:  Statistical  Aspects  and  Interaction  with  Squeezed  Light,”  Physical 
Review  A,  55:1634-1647  (February  1987). 

2.  Ackerhalt,  J.  R.  and  J.  HY.  Eberly.  “Quantum  Electrodynamics  and  Radia¬ 
tion  Reaction:  Nonrelativistic  Atomic  Frequency  Shifts  and  Lifetimes,”  Physical 
Review  D,  10{\Qi)  (November  1974). 

3.  Bagnell,  Richard  J.  Design  and  Characterization  of  Optically  Pumped  Vertical 
Cavity  Surface  Emitting  Lasers.  MS  thesis,  Air  Force  Institute  of  Technology, 
1992. 

4.  Benivenga,  G.  and  A.  Messina.  “New  Quantum  Effects  in  the  Dynamics  of  a 
Two-Mode  Field  Coupled  to  a  Two-Level  Atom,”  Journal  of  Modem  Optics, 
41  (5)  (1994). 

5.  Bjork,  Gunnar,  et  al.  “Modification  of  Spontaneous  Emission  Rate  in  Planar 
Dielectric  Microcavity  Sturctures,”  Physical  Review  A,  44  W  (July  1991). 

6.  Book,  David  L.  “NRL  Plasma  Formulary,”.  Naval  Research  Laboratory,  Wash¬ 
ington,  DC,  1990. 

7.  Bransden,  B.  H.  and  C.  J.  Joachain.  Physics  of  Atoms  and  Molecules.  New 
York:  John  Wiley  and  Sons  Inc.,  1995. 

8.  Cohen-Tannoudji,  Claude,  et  aJ.  Quantum  Mechanics,  Vol.  I.  New  York:  John 
Wiley  and  Sons,  1977. 

9.  Gou,  Shih-Chuan.  “Quantum  Behavior  of  a  Two-Level  Atom  Interacting  with 
Two  Modes  of  Light  in  a  Cavity,”  Physical  Review  A,  40(9)  (November  1989). 

10.  Guo,  Guang-Can  and  Shi-Biao  Zheng.  “Generation  of  Schrodinger  Cat  States 
via  the  Jaynes-Cummings  Model  with  Large  Detuning,”  Physics  Letters  A,  223 
(December  1996). 

11.  Hillery,  Mark  and  Robert  John  Schwartz.  “Time- Averaged  Properties  of  the 
Jaynes-Cummings  Model:  Effects  of  Detuning,”  Physical  Review  A,  43{S) 
(February  1991). 

12.  Janowicz,  Maciej  W.  and  J.  M.  Ashbourn.  “Dynamics  of  the  Four-Level  A 
System  in  a  Two-Mode  Cavity,”  Physical  Review  A,  55(3)  (March  1997). 

13.  Jaynes,  E.  T.  and  F.  W.  Cummings.  “Comparison  of  Quantum  and  Semiclassical 
Radiation  Theories  with  Application  to  the  Beam  Maser.”  Proceedings  of  the 
IEEE 51.  89-109.  January  1963. 


BIB-1 


14.  Kressel,  Henry  and  J.K.  Butler.  Semiconductor  Lasers  and  Heterojunction 
LEDs.  Quantum  Electronics  Series,  Academic  Press,  Inc.,  1977. 

15.  Loudon,  R.  The  Quantum  Theory  of  Light  (2nd  Edition).  Oxford:  Clarendon 
Press,  1983. 

16.  Marion,  Jerry  B.  and  Mark  A.  Heald.  Classical  Electromagnetic  Radiation  (2nd 
Edition).  Harcourt  Brace  Jovanovich,  1980. 

17.  Meystre,  P.  and  M.  Sargent  III.  Elements  of  Quantum  Optics  (Second  Edition). 
New  York:  Springer- Verlag,  1991. 

18.  Noble,  Michael  J.  “Physics  of  Surface  Emitting  Lasers.”  Dissertation  prospectus. 
Air  Force  Institute  of  Technology,  Wright- Patterson  AFB  OH,  September  1996. 

19.  Schlicher,  Rainer  R.  “Jaynes-Cummings  Model  with  Atomic  Motion,”  Optics 
Communications,  70{2)  (February  1989). 

20.  Shore,  Bruce  W.  and  Peter  L.  Knight.  “The  Jayne-Cummings  Model,”  Journal 
of  Modem  Optics,  .^f?:1195-1238  (1993). 

21.  Takahashi,  Ichiro  and  Kikuo  Ujihara.  “Quantum  Mechanical  Description  of 
Spontaneous  Emission  in  a  Microcavity  in  Terms  of  Admittance,”  Physica  B, 
227  (1996). 

22.  Tavis,  Michael  and  Frederick  W.  Cummings.  “Exact  Solution  for  an  N-Molecule- 
Radiation-Field  Hamiltonian,”  Physical  Review,  170:379-384  (June  1968). 

23.  Towe,  Elias.  A  Research  Report  on  Vertical-Cavity  Surface-Emitting  Lasers. 
Technical  Report  RL-TR-95-97,  Air  Force  Materiel  Command  Griffiss  Air  Force 
Base,  New  York:  Rome  Laboratory,  September  1995  (AAL-9631). 

24.  United  States  Air  Force  Scientific  Advisory  Board.  New  World  Vistas:  Air  and 
Space  Power  for  the  21st  Century:  Directed  Energy  Volume.  Washington,  D.  C.: 
USAF  Scientific  Advisory  Board,  1995  (AAM-4728). 

25.  Vurgaftman,  Igor  and  Jasprit  Singh.  “Semiconductor  Microcavity  Light- 
Emitting  Structures.”  Proceedings  of  the  SPIE2397.  515-525.  1995. 

26.  Wu,  Ying  and  Xiaoxue  Yang.  “Jaynes-Cummings  Model  for  a  Trapped  Ion  in 
Any  Position  of  a  Standing  Wave,”  Physical  Review  Letters,  78{16)  (April  1997). 

27.  Xie,  Y.  B.  “Two-Level  Atom  and  Multichromatic  Waves  in  a  Lossless  Cavity,” 
Journal  of  Modem  Optics,  44i^)  (1997). 

28.  Yamamoto,  Y.,  S.  Machida  and  G.  Bjork.  “Microcavity  Semiconductor  Laser 
with  Enhanced  Spontaneos  Emission,”  Physical  Review  A,  .^.^(1)  (July  1991). 

29.  Yang,  Xiaoxue  and  Yuanjie  Li.  “Unified  and  Standardized  Procedure  to  Solve 
Various  Nonlinear  Jaynes-Cummings  Models,”  Physical  Review  A,  55(6)  (June 
1997). 


BIB-2 


Vita 


Capt  Dustin  Ziegler  was  born  in  Burlington,  Vermont  in  1971.  He  attended 
Mount  Mansfield  Union  High  School  in  Jericho,  Vermont,  and  graduated  in  the  spring 
of  1989.  Later  that  year,  he  enrolled  at  the  University  of  Vermont  in  Burlington  on  an 
Air  Force  Reserve  Officer  Training  Corps  scholarship  through  a  cross-town  program 
with  St.  Michael’s  College  in  Colchester,  Vermont.  In  1993,  he  graduated  Cum 
Laude  with  a  Bachelor  of  Science  Degree  in  Physics  and  received  his  commission  as 
a  Second  Lieutenant. 

His  first  assignment  was  at  Los  Angeles  Air  Force  Base  where  he  worked  in 
the  Space  Test  and  Evaluation  Directorate  as  the  Satellite  Payload  Engineer  for  the 
Advanced  Research  and  Global  Observation  Satellite.  In  June  of  1996,  he  entered 
study  at  the  Air  Force  Institute  of  Technology. 


Permanent  address;  HCR  Box  227 

Richmond,  Vermont  05477 


VITA-1 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNo.  0704‘0188 


1204,  Arlington,  VA  22202-4302,  and  to  the  Cjffice  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


December  97 


4.  TITLE  AND  SUBTITLE 

SPONTANEOUS  EMISSION  IN  MICROCAVITY  LASERS 


6.  AUTHOR(S) 

Dustin  P.  Ziegler,  Capt,  USAF 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Institute  of  Technology 

2750  P  Street 

WPAFB  OH  45433-7765 


Master's  Thesis 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GAP/ENP/97D-1 1 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Dr  John  P.  Loehr 
Heterojunction  Physics  Branch 
Avionics  Directorate  (WL/AADP) 

Wright  Patterson  AFB  OH  45433-7765 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTE 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 
APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

An  understanding  of  spontaneous  emission  processes  within  microcavities  is  crucial  in  addressing  the  need  to  make  tomorrow's 
microlasers  more  efficient.  One  approach  to  improving  the  device  efficiency  is  to  reduce  the  threshold  input  energy  at  which 
lasing  begins  to  occur.  It  has  been  suggested  that  the  threshold  in  a  microcavity  laser  can  be  decreased  by  increasing  the 
fraction  of  spontaneous  emission  into  the  lasing  mode,  this  can  be  accomplished  by  preferentially  coupling  the  gain  medium  of 
the  laser  to  the  electromagnetic  cavity  mode  of  interest.  It  therefore  becomes  necessary  to  understand  the  mechanism  by 
which  this  coupling  takes  place.  This  research  develops  a  fully  quantum  mechanical  description  of  the  interaction  between  a 
gain  medium  modeled  as  a  two  level  atom  and  a  multimode  electromagnetic  field  in  a  microcavity.  Atomic  transitioii 
probabilities  are  computed  for  systems  in  which  the  atom  couples  through  a  single  photon  process  to  electromagnetic  cavity 
modes  which  range  in  number  from  two  to  2000.  Calculations  performed  for  cavities  with  widely  spaced  modes  demonstrate 
that  atoms  exhibit  Jaynes-Cummings  behavior  when  closely  tuned  to  one  mode.  Detuning  of  the  atom  from  the  mode  inhibits 
the  exchange  of  energy,  while  increasing  the  strength  of  the  coupling  to  the  mode  amplifies  this  exchange.  Two  level  systems 
strongly  coupled  to  many  closely  spaced  modes  exhibit  spontaneous  emission  rates  characteristic  of  an  atom  in  free  space.  ‘ 


14.  SUBJECT  TERMS  15.  NUMBER  OF  PAGES^ 

Quantum  Electrodynamics,  Field  Theory,  Diode  Lasers,  Laser  Cavities,  Spontaneous  Emission  1 19 

16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICAtlON  bo.  LIMITATION  OF  ABSTRAC 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


Unclassified 


Unclassified 


Unclassified 


Designed  using  Perform  Pro,  WHS/DIO=),  Oct  94 


